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Abstract 

Motivated by recent studies by Dorey, Pocklington and TateojSH ESjfor uni- 
tary minimal models perturbed by (/'i,2 5 we examine the thermodynamics of one 
dimensional quantum systems, whose counterparts in the 2D classical model are 
the dilute Al models in regime 2. The functional relations for arbitrary values of L 
are established. Guided by numerical evidences, we obtain a set of coupled integral 
equations from the established relations, which yields the evaluation of the free en- 
ergy at arbitrary temperature. In the scaling limit, the integral equations coincide 
with the thermodynamic Bethe ansatz equations (TBA) proposed in thereby 
support their results. The new Fermionic representations of the Virasoro characters 
are shortly remarked. 
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1 Introduction 



The Andrews-Baxter Forrester model T] has been a prototype in the studies of integrable 
lattice models as well as the integrable $1,3 perturbation of the minimal unitary CFT. 

The dilute model proposed in [3 EI deserves the same attention, as it represents a 
lattice analogue of the 0i 2 and 02,i perturbations of the minimal unitary CFT. 

Especially, there are abundant research results on the dilute ^3 model in regime 2, 
which corresponds to the model of high interest, M3^4 + $1^2 theory. The latter belongs to 
the same universality class as the Ising model in a magnetic field at Tc, and is characterized 
by the famous Eg structure|3]- J2]- There are several evidences for this coincidence, the 
central charge [21 E], surface exponents [12] universal amplitude ratio O [13 HEl scaling 
dimensions jl7j and excitation spectra ^Hl- The most impressive demonstration may be 
presented in ^H], which confirms the hidden Eg structure behind the thermodynamic 
Bethe ansatz equation (TBA, for short) of the dilute ^3 model. See also [213 [Zlj- 

Similarly, M4^5-|-$i^2(^6,7+'^'i,2 ) theory has the hidden Ey{Eq) structure behind In 
spite of lack of relevant string hypothesis, the technique of the quantum transfer matrix 
(QTM) and the fusion relations make it possible to manifest the desired Lie algebraic 
structure in the TBA of the dilute ^3,4,6 in regime 2[22j l23]. 

The aim of the present report is to extend the result further; we will establish a closed 
set of fusion relations of the dilute Al model for arbitrary L. From these relations, we 
derive the nonlinear integral equations (TBA) which yields the evaluation of the free en- 
ergy of the corresponding ID quantum model at arbitrary temperatures. This is obviously 
motivated by the recent progress in the study of the models as perturbed conformal field 
theories [211 122] • The perturbation theory by an operator $12 or $2,1 encountered serious 
technical difficulties. The systematic studies on the bootstrap procedure on S matrix have 
been initiated in [22] and |2Z|- See also [23 1211 • The latter approach bases on the scaling 
q— state Potts field theory and will be relevant in our view point. It is only recently that 
the closed bootstrap procedure, associated to the Potts field theory, has been obtained for 
a set of S— matrices for a wide range of parameters by Dorey, Pocklington and Tateo^lj. 
The check of the results against a finite size system however, again poses a challenging 
problem because of the non-diagonal nature of the scattering theory. In |25], a set of in- 
volved TBA is conjectured by ingenious insights from special cases which have similarity 
to the TBA for the sine-Gordon model. They prepare pieces of TBA equations, one of 
which comes from Zn models[3Q], then glue them together by respecting symmetries and 
by physical requirements. The resultant equations pass several nontrivial checks, however, 
lack firm grounds. It is thus interesting to see if it is possible to recover the Y— system 
[SH IS21 Hn] ; the algebraic version of the TBA, from the fusion relations of the dilute A^ 
model. 

The functional relations among transfer matrices (the T— system) are determined 
by the fusion relations [SH ISE]- The decomposition of the product of two modules at the 
singular points of the Boltzmann weight rules the latter. Thus the problem seems to 
be straightforward, at the first sight. This is , unfortunately, not true. There are often 
infinitely many functional relations for a given model. One thus has to choose a subset 
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which offers the closed T— system with finitely many entries. Moreover, among the sets, 
the only those of which both sides contain fusion transfer matrices with good analytic 
properties (to be specified) will be relevant for the TBA. Therefore the problem of finding 
the relevant functional relations so as to reach the TBA is far from trivial. 

It is well known that the dilute Al model has Uq{A2 ) symmetry at the critical point. 
The off-critical extension, deformed Virasoro algebra Vx,r{^^2^) plays a role in the vertex 
operator approach jHZj. This elegant result, however, is not of our direct relevance. The 
crucial fact for the present study is, rather, that the fusion structure which is valid even 
away from criticality [SB]; it possesses the fusion structures similar to Uq{A^^'') and Uq{A^2^) 
. We will relate the former to breathers while the latter to the pseudo-particles constituting 
ladders in their TBA diagrams. In addition, we will introduce a set of objects which we call 
"magnon-like" , characterized by similar diagrams to those associated with simply- laced 
Lie algebra (A, D and E). It turns out that they constitute closed functional relations with 
the good analyticity. They are three almost independent relations but for one common 
object Teg, which plays the role of glue. Thus the construction employed here makes the 
meaning of gluing in very explicit. 

We will introduce auxiliary functions, fusion transfer matrices, ^ labelled by skew 
Young diagrams|Sni HHl 1^ In view of the latter diagrams, proofs of the functional 
relations seem to be similar to those of the decomposition rules in the symmetric group. 
There are, however, extra complication due to the fact that our coupling constant in a 
sense lies at a "root of unity". Sometimes two (or even more) diagrams, with completely 
different disguise, correspond to an identical fusion transfer matrix. One must choose 
relevant expressions (diagrams) case by case in order to obtain suitable decompositions. 
This makes the proof of the functional relations simple but lengthy, unfortunately. There 
are four families of functional relations depending on the values of L and they need 
separate treatment at the present stage. 

The validity of the obtained algebraic equations does not necessary legitimate the 
integral equations (TBA). We need another step to clarify the analytic structure of Y— 
functions. This remains an important issue to be understood. In this report, we employ 
the direct numerical calculation by various range of parameters and adopt a simple con- 
jecture on the analyticity. With the established functional relations and the assumption 
on the analyticity, it is then straightforward to transform the T— system into Y— system 
and then Y— system into coupled integral equations S35^. In the scaling limit, the TBA 
proposed in is recovered from them. A simple ansatz, suggested by numerics, on the 
zeros of fusion transfer matrices remarkably leads to the correct mass ratio of the particles 
determined by S— matrix poles. 

This paper is organized as follows. In the next section, we give a brief review on the 
dilute Al models and the QTM method. Some obvious consistency of the lattice models 
and the CFT perturbed by $1,2 will be remarked. Section 3 is devoted to the discussion on 
fusion transfer matrices which are parameterized by skew Young diagrams. We introduce 

^We expect them to be the eigenvalues of transfer matrices of which fusion types are described by the 
diagrams. This is however yet to be proved. As our argument does not rely on whether this coincidence 
is true or not, we simply call the auxiliary functions as (eigenvalues of) fusion transfer matrices. 
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there important kinds of objects in section 4, 5 and 6. The first object introduced in 
section 4 seems to correspond to the breathers in [21] while the second one in section 5 
is responsible for the ladder structure in the Y— system diagram in jSHI- The section 6 is 
divided into 4 subsections, each is devoted to the discussion on "magnon-like" objects. In 
case of the dilute Al model, L even, a fundamental role seems to be played by a "kink " 
transfer matrix, which is introduced in section 7. We will give a remark on the "reality 
property" of the magnon-like objects in section 8. Once fusion relations, the T— system, 
are established, it is readily transformed into Y— system. The explicit procedure will be 
discussed in section 9. The proofs of the T— system (or the "magnon-like" t— system) are 
simple but lengthy. Those for L even (odd) are somewhat similar. For brevity, we thus 
present the proof for L = 4/c + 2 in the section 10 and L = 4A; -|- 1 in the appendix F. In 
section 11, it will be shown that the conjectured TBA is naturally restored in a scaling 
limit. We conclude the paper with brief summary and discussion in section 12. A new 
fermionic formula for the Virasoro character found for M^^q is briefly commented. 

2 The dilute Ai models and the quantum transfer 
matrix 

The dilute Al model is proposed in [21 as an elliptic extension of the Izergin-Korepin model 
The model is of the restricted SOS type with local variables G {1,2, ■■ ■ ,L}. The 
variables {a,b} on neighboring sites should satisfy the adjacency condition, |a — 6| < 1, 
which is often described by a graph in figO The following RSOS weights are found in j2] 




Figure 1: An incidence diagram for the dilute A^ model. The local states corresponding 
to connected nodes can be located to nearest neighbor sites on a square lattice. 

which satisfy the Yang-Baxter relation. They are parameterized by the spectral parameter 
u, the elliptic nome q = e~^. 
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00 

A{x,r) = 2gi/Sma;JJ(l-2g2"cos2a; + g^'^)(l-g2'^), 
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00 
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q - ( ^^a^J!M. 

- ^ ^> e,{2ay 

Note an important property, the standard initial condition, that only those weights with 
identical variables on the NE and the SW diagonal remain nonzero when m — 0. 

The crossing parameter A needs to be a function of L for the restriction. The model 
exhibits four different physical regimes depending on parameters j21 Ej, 

• regime 1. < m < 3, A = 4(f^, L>2 

• regime 2. < m < 3, A = -^^ > 3 



• regime 3. 3-f<M<0, A = L > 3 

• regime 4. 3-f<M<0, A = ^^f^, L>2. 

In the low temperature limit, the ground states correspond to ferromagnetic ordered states 
in regime 1 and 2 in which we are interested. 

As usual, we define the row to row transfer matrix Tj^tr(^) by 

M 

Denoting its largest eigenvalue by the same symbol, Tj^trI""); the free energy per site of 
the 2D classical model in thermodynamic limit reads, Pf = — limM^oo ]g log^RTR(^)- 
Several other quantities other than the free energy have been evaluated for the models. 
For example, we mention the explicit one point functions elaborated in [T7] . 

We are not dealing with the dilute Al model directly but its one dimensional quantum 
counterpart. 

The Hamiltonian for the latter is defined by 

H = -2|-lnTRTR(«)L=o (2) 

as in pni- The explicit form of the Hamiltonian is not needed in the following argument. 

We are interested in the free energy of the ID model at finite temperatures. The 
success in the previous studies on the thermodynamics of ID quantum systems depends 
on the efficient conjectures on the dominant solutions of the Bethe ansatz equation in 
the thermodynamic limit , the string hypothesis. It seems too difficult to perform the 
same procedure for arbitrary L. Instead, we apply another machinery, the method of the 
quantum transfer matrix (QTM)|^-|48j. 

The quantity of our interest is the partition function of the ID model. 



^ID 



tre"^^. 



The first step is to rewrite this as the partition function of the 2D model, defined on 
M X N square lattice by a simple trick, 

= lim tr(e"^^)^ = lim Tr(r^)^ = lim TrZao- (3) 

N^oo N^oo N—*oo 

The fictitious dimension, N, is referred to as the Trotter number. is a kind of transfer 
matrix , composed of matrix elements of e^iv^. 

In order to relate to the row to row transfer matrix above, we introduce a rotated 
transfer matrix, Tj^trI"") > 

M 



1 (^j+i 
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Thank to ^ and the standard initial condition, we immediately see, 



where Pr{Pl) is the right (left) shift operator, commuting with the Hamiltonian and 
PlPr = PrPl = the identity operator. Then the desired expression for is obtained 

= ^rtr(^)^rtr(^) 

where the limit ^ oo is assumed in the both sides. The eigenvalues of are highly 
degenerated in ^ oo , which makes the summation in the partition function difficult to 
perform. The intriguing observation ^3] is that if one defines yet another transfer matrix 
T, the quantum transfer matrix (QTM), propagating in the horizontal direction, then 
there is a gap in between the largest and the second largest eigenvalues which remains 
non vanishing in the limit N ^ oo. Now that Zid ~ TrT*^ and M — >■ cxd, the only largest 
eigenvalue of T is necessary in the evaluation of (jHI). 

In the most sophisticated formulation |49 j . the QTM is defined with one further variable 

X, 



(7QTM(«,^)){a} = n 

By taking a; = one recovers T . 

The commutative property of QTMs with fixed u is crucial in the QTM formulation, 

[^QTM(w,a;),TQTMKa;')] = 0. 

It represents the existence of infinitely many conserved quantities commuting with QTM, 
which may be a key in the exact enumeration of the free energy at arbitrary temperatures. 
Since we will adopt identical value for m = for commuting QTM, it will be dropped 
afterwards. 

Let us explicitly write again that the free energy per site is represented only by the 
largest eigenvalue of Tq^pj^ with x = and u = 

Pf = — lim — InTr exp(—d7i) 





b2j 


u+ix 






127 





b2j 


T 




02, 



N~ 



lim In (the largest eigenvalue of TqtmI^^ = 0) 



Our aim is thus to evaluate the largest eigenvalue of QTM. This is again not so simple 
problem as it seems at the first sight. The coupling constant u depends on the fictitious 
system size. Trotter number, which makes it difficult to take the limit oo. 

Our strategy is to introduce auxiliary transfer matrices (fusion transfer matrices) com- 
muting with QTM and to explore functional relations among them j22l 123 EHl IHH IHH tHB] • 
We find the latter relations, if suitable subset is chosen, are efficient enough to yields the 
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evaluation of the largest eigenvalue of Tq^^i for arbitrary N: The functional relations can 
be transformed into coupled integral equations whose solution yields the desired quantity. 
The resultant equations admit to take the limit N —>■ oo without any difficulty. 

In the next section, we will introduce the fusion transfer matrices. Before doing so, 
we mention a simple evidence that the dilute Al model (or its ID analogue) is a lattice 
analogue of the minimal unitary model perturbed by $1^2 or $2,1 operator. In TBA for- 
mulation, the dressed energy function e and Y = exp(e) are main objects. It is pointed 
out in [23] that, within framework of TBA, it seems too difficult to show analytically that 
the small r expansion coincides with the result from the conformal perturbation theory. 
The hidden periodicity of the Y— system, found by a brute-force substitution, partially 
resolves the problem ^T] . 

We will show, within our formalism, Y is expressible by ratios of fusion transfer ma- 
trices. Therefore in our view point, the periodicity of Y (and e ) is a consequence of the 
periodicity of the original Boltzmann weights. In view of the rapidity variable 6 = the 
periods read 

_ . ll^gll^, for regime 1 

- S 4(L+l)7r r 



3{L+2) ' 



for regime 2. 



From this one can read off the scaling dimension of the perturbed field by A = 1 — p-, 
explicitly. 



p+3 



for regime 1 



A = <^ ^ (5) 

lifei)' for regime 2. 
Note that we adopted a different identification, 

{L + 1 for regime 1 
T f ■ O 

L for regime 2. 

They coincide with the scaling dimensions A2.1 and Ai 2 in the unitary model Mpp+i. 
This thus supports the hypothesis that the dilute Al models in regime 1 (regime 2) are 
lattice regularizations of perturbed minimal models Mp^p+i under the identification (jHl). 

From now on, we restrict our argument to the regime 2. As is announced, the fusion 
transfer matrices will be examined in the next section. 



3 The fusion transfer matrices and quantum Jacobi- 
Trudi formula 

Any model integrable in the sense of the Yang-Baxter relation offers a systematic way of 
generalizations, the fusion procedure. It utilizes the singular values of spectral parameter 
at which the Boltzmann weights reduce to projectors to a subspace. In IHHl, a detailed 
description is given for the explicit procedures to obtain "generalized" Boltzmann weights. 
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The explicit construction of these weights for various fusion types is not of issue in this 
report. Rather, the main objects of our interest are functional relations which involves 
QTM. In the following we introduce several auxiliary functions which generalizes QTM. 
We expect that they corresponds to eigenvalues of transfer matrices made by the fusion 
procedure. We, however, do not attempt to prove this, as it is not necessary for the 
evaluation of free energy. 

The most important building block in our argument is the following expression for the 
(any) eigenvalues of QTM, which we will denote by Ti{x), 

A.( ^3..,. , \ .^ Q{x-\i) , ^, 3 3 Q(x-fz)Q(x + |z) 



+ w <P[x — -t)(p{x — -I)— j— -, (7) 



N 

Q{x) := JJ^ h[x — Xj] 

A( \ ( h[x + {l-u)i]h[x-{l-u)i] ^^ 

where w = exp{i^^) and i = 1 for the largest eigenvalue sector. It is often referred to as 
the dressed vacuum form. 

The parameters, {xj} are solutions to "Bethe ansatz equation" (BAE)[Tn]. 

^ Qjxj - i)Q{xj + 2i) ^. _ ^ . . . ^ 

(f){xj - i) Q{xj + i)Q{xj - 2i) ' 

We introduce 

j^iL+ll^ for regime 1 

^ - i 4(L+1) r • o 

I ' , for regime /. 

Then the periodicity Ti{x) = Ti[x + Pi) follows obviously from the above expression. 
We represent the three terms in eigenvalue of the quantum transfer matrix by three 
boxes with letter 1,2 and 3 suffixed by the spectral parameter . 



2'^^ 2'Q{x-\{) 

^/ 3.,^, 3.,Q(x-|z)Q(x + f^) 
2 := d)ix+-i)d)ix 1)—, H — rr (10) 

\^ = w 0(a; - -i)(t){x - -i)- 



2'^' 2'Q{x + Uy 
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These boxes are analogue to those for Young tableaux and play a fundamental role in 
the following. 



3.1 The symmetric fusion structure 



In the bootstrap procedure of the S— matrix, (p^ property is most essential. The similar 

(2) 

feature , referred to as the 05 fusion structure, as been discussed in the dilute Al model 
|38j . This comes from the singularity of the RSOS weights at u = ±2. 

Due to the fusion procedure, the eigenvalue of a fusion QTM can be represented by 
sum over products of "boxes" with different letters and spectral parameters. The point 
is that the sum of such boxes can be identified with Semi-Standard Young Tableaux 
(SST) for s/3. Let us present a simple example. See [SHI Ell EH ESI ES] for details. By the 
fusion procedure, one can construct a transfer matrix of which auxiliary space acts on a 
symmetric subspace of V x V. We associate this to the set of the SST , iii2, {ii < «2)- 
The eigenvalue of the transfer matrix is then represented by. 



x+i 



il<i2 



In a same manner, one can introduce fusion models based on general Young diagrams, of 
which eigenvalues are specified by their shapes. On each diagram, the spectral parameter 
changes by +2i from left to right and —2i from top to the bottom. 

First, we argue QTMs associated to rectangular shape diagrams. Although they are 
simpler, some of the properties below are crucial for the discussion for QTMs associated 
to complex Young diagrams. First, due to identities. 



X + l 

x—i 



x+i 
x—i 



02(xM, 



x+t 
x—i 



02 



x+2i 



x-2i 



(f)2ix) := (j){x ± -i) 



(11) 

(12) 



the QTMs from 2 x m ( 3 x m) Young diagram can be reduced to those from 1 x m (or 
just scalars). Hereafter we adopt abbreviations for any function /, 

/(x ± ia) := f{x + ia)f{x — ia). 

(2) 

The first three relations brings additional symmetry besides s/3, and make the 03 

(2) 

symmetry explicit. We call this the 03 structure hereafter. 

Second, the eigenvalues of 1 x m fusion QTMs have the "duality" in the following 
sense. Let us denote a renormalized 1 x m fusion QTMs by Tm(x); 



Trr,(x) 



fn 



(X 



El 





i2 







(13) 
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where the semi-standard (SS) condition < ^2 < • ■ ■ < "^m is imposed on the summation. 
The spectral parameters are assigned x — i{m — 1) ■ ■ ■ x + i{m — 1) from left to right. The 
renormalization factor, which is the common factor of the expressions from tableaux of 
length m, is given by 

"i^ 2m - 1 

such that the resultant T^'s are all degree 2N w.r.t. h[x + shift ]. For later convenience, 
we also denote by T^°'^(x), the un-renormalized Tm{x), namely, the object without the 
factor fm{x) in (fT^ . Obviously, we have a periodicity due to Boltzmann weights; 

Tm{x + Pi) = Tm{x). (Tr"(x + Pi) = rr"(x)) 

where P is defined in ©. 

Hereafter, we frequently use functions of which argument is shifted by the half period in 
the imaginary direction, i.e., x + i^. We shall denote them by a symbol ^. For example. 



p 

T^ix) =Tmix + —i). 



2 

(2) 

The following functional relations are direct consequence of the structure. 

Tm{x - i)Tm{x + i) = gm{x)Tm{x) +Tm+i{x)Tm-i{x), m > 1 (14) 

3 

g-mix) = (j){x±i{m + -)), 
T_i(x) := 

To{x) := f2{x). (15) 



The scalar function gm{x) satisfies, 



g2L-l-m[X) 



gm{x), if L is even 
g^ix), if L is odd. 



The adjacency condition leads to T2l{x) = T2l+i{x) = 0. Thus one deduces the duality 
relations, 

y . s _ \T2L-i-ni{x), m = 0, ■ ■ ■ , 2L for L even 
" l7^2Vl-„^(a;), ^ = 0,--- ,2L for L odd 

for the solutions to eq. lfT^ . Note that this duality is not valid for T°°'^(x); obviously 
T2L-i-mi^) '^m^i"^) havc different order in h[x]. These relations can be in principle 
proved by using explicit fusion weights. We have at least checked this numerically and 
assume the validity in this report. 

The duality brings about many relations. It thus happens seemingly quite different 
equations are identical; an identical object can have (infinitely ) many equivalent expres- 
sions. It sometimes makes the proof of the functional relations quite involved. 
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Before closing this section, we make a comment. 

The simple relation (fT^ suggests a choice of Y function, 



Tm+i{x)Tm-i{x) 

gm{x)Tm{x) 



(m = I,-- - ,2L-2). 



(17) 



This leads to a set of simple equations (the Y— system) 



ym{.X - i)ym{.X + i) 



(i + y^+i(x))(i + y^-i(x)) 

1 + ym{x) 



(18) 



The relations seem to be tractable and they indeed work well in a nonunitary case 
|57j . This is not the case for the unitary models. The transformation of the Y— system 
into a simple set of integral equations requires the both sides to be Analytic and NonZero 
and to have Constant asymptotic behavior as |a;| ^ oo (ANZC) in certain strips. 

As is demonstrated numerically in appendix in the case of L = 8 in appendix [Gl this 
requirement is not satisfied for (fTHj) . Therefore we must seek for alternatives, which may 
not be found only among symmetric fusion transfer matrices. 

In the next subsection, we thus consider a much wider class of skew Young diagrams. 

3.2 The quantum Jacobi-Trudi formula 

Let /i and A be a pair of Young diagrams satisfying /ij > Xi,Wi. We subtract a diagram A 
from /i. The resultant "narrower" one, consisted of (yUi — Ai, /i2 — A2, ■ ■ ■ ) boxes is called a 
skew Young diagram fi — X. (The usual Young diagram is the special case that A is empty, 
and we will omit A in the case hereafter.) 

Consider a set of semi-standard skew Young table of the shape fi — X. The spectral 
parameters are assigned to boxes in a table such that they change by +2i from left to 
right and —2i from top to the bottom. The spectral parameters of the head (left-top) and 
the tail (right-bottom) are fixed to be x + i(/ii + fi[ — 2) and x — i{fii + fi[ — 2) where fi'i 
denotes the depth of the diagram. One identifies each box in a table with an expression 
under the rule ()1U|) with a shift of the spectral parameter. Then the product over all 
constituting boxes yields the desired expression for a table. We then have the following 
theorem which is a quantum analogue to Jacobi-Trudi formula for the Schur function. 

Theorem 1. Let T^/x{x) be the sum over the resultant expressions divided by a common 
factor, Y[%i ffij-Xjix+i{fi[—fii+^j+Xj—2j+l)). Then the following equality holds lS^ IJUij. 

Tfji/x{x) = det i<j^k<i,[iTf,^~Xk-j+kix + i{fi[ ~ fii + fXj + Xk - j - k + 1))) (19) 
where Tm<o := 0. 

T^/x{x) may be naturally identified with the eigenvalue of QTM corresponding to 
fusion /i — A, although this is yet to be proved. We note that the pole-free property of 
T^/x{x), obvious from the quantum Jacobi-Trudi formula. 
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Due to the semi-standard condition, products of T^/x{x) enjoy the decomposition sim- 
ilar to that of the products of Young diagrams. This observation is crucial in the proving 
the functional relations below. One, should , however note that each box can only stay at 
positions with identical spectral parameter. Therefore, the product of two T shows non- 
trivial decomposition only when corresponding diagrams are situated at proper relative 
positions specified by their spectral parameters. 

For illustration, let us prove (fT^ diagrammatically. We consider T^'^{x + i)T^^{x — i). 
The spectral parameters are assigned for each box, x — {m — 2)i, ■ ■ ■ ,x + mi from the 
left to the right for T^'^{x + i) , and x — mi, ■ ■ ■ ,x + {m — 2)i ioi T^"^{x — i). Thus two 
diagrams are situated so that these two width m rectangles piled vertically (fig. |21). @ 

We consider a particular assignment of letters, ii, ■ ■ ■ ,im for the upper and ji, ■ ■ ■ ,jm 
for the lower diagram. T^°'^(x + i)T^'^{x — i) is obtained by the sum over tableaux under 
the condition, ii < i2 < ■ ■ ■ < im and ji < j2 < ■ • • < jm, respectively. 



ii 


^2 


13 




^5 
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17 


is 


^1 




^3 


^4 


^5 




^7 


^8 



Figure 2: A typical table in the product T^^{x+i)T^^{x—i). Hereafter, disjoint diagrams 
represent their product. 

The first term in the Ihs of (|T^ comes from the configurations, ik < jk, {k = 1, ■ ■ ■ ,m). 
The resultant table is of the shape 2 x m and satisfy the semi-standard condition. Due 
to the reduction (fTTj). it is a member of T^°"(x) times HJLi'i'^zla^ — + 1) + 2ij). 
The second term comes from the configurations which breaks the condition (fig. By 
ia > ja, we mean the first pairs breaking the condition from the left end . Then the 
ordering within z's and j's concludes that the new sequences ji, ■ ■ ■ , ja, ia, ia+i, ■ ■ ■ ,im and 
ii, - ■ ■ ,ia-i,ja+i, ■ ■ ■ ,jm satisfics the semi-standard condition. One can also check that 
the new sequences satisfy the proper spectral parameter assignment, namely increment 
of the spectral parameter by 2i from the left to right. Thus after the recombination of 
sequences , we obtain T^_^i{x)T^_^^{x) . Finally, by dividing the both sides by the factor 
fm{x ± i) , one arrives at (fT^ . 

The above procedure is formally interpreted as follows. Place two diagrams properly 
according to their spectral parameters. We just join them or recombine them such that 
the total set of spectral parameters is conserved. 

The case studies for L = 3,4, 6, however, suggest that T^/\{x) is not the right object 
in the discussion of functional relations. Instead, we should introduce K^i\{x), which is 
analytic under BAE, from T^jx^x) by putting Tm>2L{x) = in the latter, 

A^/a(x) := T^/xix)/.{T^>2L{x) = 0}. (20) 

The pole-free property of A^/x{x) is obvious from (fT^. 

The diagrammatic decomposition applies to T^/x{x) but not to A^/a(x) in a strict 
sense. Therefore the above decomposition rules, demonstrated above, may not seem to 
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Figure 3: The pair (i4,j4) is the first from the left which breaks the semi- standard 
condition. Then the recombination of tableaux occurs there, into hatched and unhatched 
ones. 



be helpful in our discussion. For most of case which we treat below, however, the skew 
diagram — A is no so huge and A and T coincide. Thus the above rules are quite often 
applicable and play significant roles ; as already remarked, there are plenty of expressions 
for an identical object. Therefore it is sometimes crucial to choose "the most" relevant 
expression among many. The above graphical decomposition provides hints on the optimal 
candidate in proving the functional relations . 

Before closing the section, we mention an important diagrammatic symmetry of A^/x 



Corollary 1. Let Y and Y be skew diagrams. Y is obtained from Y after 18IT rotation. 
The following relation is a consequence of the eq M9^] 

Ty{x) = {Ty{x*)y (21) 

where the asterisk stands for complex conjugations. 

The proof is direct; one substitutes Kj = (/ii — \d+i-i) — {j^i — jdd+i-i) for F = yU — A 
and = d. 

This corollary is crucial in establishing a property of the transfer matrices (see section 
IH}. It will be also used implicitly in the proof of functional relations in the following. 

After these preparations, we introduce the first of our three ingredients in the closed 
functional relations, the Breather T— system in the next section. 



4 Breather T— system 

The RSOS weights possesses a singularity at m = ±3 besides the u = ±2 singularity 
related to a'^^ structure jSH]- It is shown that the singularity at u = ±3 signifies the 
"hidden" s/2 structure behind the model. The breather type relation originates from this 
s/2 structure. 
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We define the "breather" transfer matrices for the dilute model in regime 2. 

Tb^(x) = Ti{x) 

Tb (x) = /TO'^2L-2,i(a^ - ^oi L odd 

\w{x)^2L-2,i{x - for L even 

Ts^ix) = ^^^^;^^A(4L-5,2L-2,2L-2)/(2L-3,2L-3)(a;) 

One immediately verifies the validity of the following functional relations, 

TBAx±i ,^~^J = To(x±i ^ ^ ^ ^ ) + (f)'' (x)TbJx) (22) 
2(L + 2)^ ^ 2(L + 2)^ V y i^3V 7 v ; 

L — 2 4 

^i?3(3^±^ 2(L + 2) ^ = To(x)To(x±i ^^^^^ ) +TB,(a:)TB,(3:). (23) 

Proof. Proof utilizes the duality (fT^ . For ()22|) , we first consider the multiplication of 
Ti(x)(= Tsj) by its dual T2i,_2(a; + (2L — l)i) . See fig. HI for a case L = 5. 
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Figure 4: The multiplication of Ti(x)(= TbJ by T2L_2(a; + (2L — The letter a in a 
box implies its spectral parameter x + ia. 

The spectral parameters are assigned so that they align horizontally. For L even (odd) 
T2L-2(a; + (2L - is equal to Tb,{x + (2L - ( Tb,{x + {2L - l)i + f i) ). Thus the 
product is equivalent to that of two Tb^. It is decomposed into two terms; the first is T2L-1 
and the second is T^g (with some shifts in spectral parameters and with coefficients.) By 
the duality, T2L-.1 is replaced by Tq. Then one arrives after the shift in the spectral 
parameter and taking account of proper coefficients. 

The second relation is shown similarly. One only should notice that dual of the diagram 
(2L — 1, 1) is given by a skew diagram(2L — 1, 2L — 1) — (2L — 2). The case L = 5 is 
depicted in fig. El . 



Figure 5: 
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Despite its disguise, the first diagram in the rhs corresponds to the products of Tq due 
to ()19|) . Then it is straightforward to verify (j^^ . □ 



By noticing 3= 



(L+2) 1"^^*^ L+2 



we understand the above relation stems from 



the singularity of the RSOS weight at m = ±3 related to s/2 structure. 

The following representation of Tb^{x) and Tb^{x) will be useful in the following ar- 
gument for the regime 2 



2(L + 2) ^ 

Tb.X^) = A(3,2)/(i)(a;). 
Proof. The first equation utilizes the representation of Tg^ (x) by the height 2 diagram. 



(24) 
(25) 



TbA^jTbA^ + 3«) = — ]-73 — A(2,i)(a; + «) + 



^>2 a; 



^2 X 



(26) 



See dnj and ^ for 02,3- 



1 


3 


-1 





Figure 6: A graphical representation of (j2SI)- The letter a in a box implies its spectral 
parameter x + ia 



We present graphically eq. ()26|l in fig. IHl By noticing 3= 



{^(mod ), and 



comparing the above with ()22|1 (after trivial shift in x), one concludes ()24j) . 
To prove the second, we consider '7^22-5 2L-2 2L-2)/(2L-3 2L-3) (^) instead of 
A(4L-5,2L-2,2L-2)/(2L-3,2L-3)(a;) in the definition of T^g. Note their difference, ((201). 
From theorem 1, '7^4i"_5 2L-2 2L-2)/(2L-3 2L-3) (^) decomposed into 3 pieces as in the case 
of usual Young diagrams. 

"^I'l^x ± 2«L)03(x) = T2ix ± 2Z)03(X)/2L-3(X ± 2iL) 

where the duality relation (fT^ is used in the second equality. 03 (x) is defined in (fT^ . 
After proper normalizations, it reads. 



^4L-5,2L-2,2L-2)/(2L-3,2L~3) (^) 



(4L-5,2L-2,2L-2)/(2L-3,2L-3) 



X 



On the other hand, the formula leads to 



= T2(x±2z). (27) 

^4L-5,2L-2,2L-2)/(2L-3,2L-3)(a^) = A(4L_5^2L-2,2L-2)/(2L-3,2L-3) (a^) +^0(2^ ±0^4L-3 (a^) • (28) 

Note T4i_3(x) = T4(x) , which is derived in parallel to the duality. Then the comparison 
(ITfll and (j2Hl) leads to 



A(4L-^5,2L-2,2L~2)/(2L-3,2L-3)(a;) = T2{x ± 2i) - To{x ± i)T4(x) = A 



(3,2)/(l) 



where fll9|l is applied in the last equality. In this manner, one verifies (OKI ] 



□ 
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The transformation of the T system to the "breather" Y— system is standard jH3 EHl 
. We define, 



0^(x)Tb3(x) 

To(x ± «2{L+2)) 

Tb,{x)Tb,{x) 
Tn(x)T()(x ± i 



(L+2) . 



Then the following "breather" Y— system is directly derived with the help of the 
"breather " functional relation 



YbA^ ± ^ 2(1+% ^ = + ^^^^^"^^^ ^^^^ 

The main goal in this report is to evaluate the largest eigenvalue of Ti{u, 0) or Tbi(0) 
in the Trotter limit. For this, our strategy is to utilize a set of closed functional relations, 
including Tb^{x). The s/2 originated functional relations are indeed enough to obtain the 
desired relations in the case of M3^4 + 0i,2- This is no longer valid for general case. In the 
next sections, we introduce further objects and relations to achieve the goal. 



5 The s/3 type functional relations 

As is mentioned in the preceding section, there is another structure analogous to the s/3 
models. Let Wm \ (a = 1, 2. m G Z>q) be a Yangian module corresponding to mA^ as a 
classical s/3 module. We associate a transfer matrix Tm\x) to win' . Then the following 
functional relations hold, 

where 7 is a parameter related to a singularity in R matrix. Tq"'' , Tm ■* and Tm ■* are given 
by trivial scalar factors times an identity operator. Thus when multiplied on a common 
eigenspace of the commuting transfer matrices, they amount to known trivial scalar func- 
tions. The above relation can be graphically represented as follows. Consider a product 
of two rectangles of the shape a x m. It decomposed into two subsets. The one is the 
products of two rectangles of the shape (a + 1) x m and the shape (a — 1) x m and the 
other is the products of rectangles , a x (m + 1) and a x (m — 1) . 

We turn back to the dilute models in regime 2. The above relations among rect- 
angular can not hold ; due to the 03 symmetry, the 2 x m rectangle is always reduced 
to the rectangle of 1 x m . Nevertheless, an analogue can be still found . Let us denote 

Td,{x) = A(j+ij)/(i)(x). 

In this notation, which is adopted hereafter, Tb^X^) in the preceding section is identical 
to Td^{x). Note also that Td^(x) = Ti{x ± 2i) and Td^^x) = To{x ± 2i). 
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Lemma 1. The following functional relations, similar to s/3 type, hold. 

T2j{x±2i) = T2(j+i){x)T2(j-i){x) +Td2,{x) (30) 

T (x±2t) = 1^^20+1) (^)^i32o-i)(^)+ ^2, (x)TB3(a;±i^y^) j even 

\TD,u.^^^)Tn,,.,,{x)+T2,{x)^^^^^^ j odd. ' 

The similarity is clearer by writing Tj^\x) = T2j{x) and Tj'^\x) = Td2^(x). The 

difference lies in that fact that (x), which amounts to a known scalar for the s/3 case, 
remains non trivial; it is given by a product of Tsg. 

Proof. The first equality (jHUj) is trivial shown by the quantum Jacobi-Trudi formula. The 
second one needs one step further. First, one uses the following relation, proved by the 
quantum Jacobi-Trudi formula, 

Td2, {x ± 2i) - Td^^^.^,) (x)Ti52U-l) i^) = T2j{x)r^2j+2,2j+l,2j)/{2,l)ix). (32) 

From the tableaux rule, a hatched part of the diagram in T(^2j+2,2j+i,2j)/{2,i){x) reduces to 
product of scalars, see fig. [3 



Figure 7: 7^2i+2,2i+i,2i)/(2,i)(a;) for j = 3. 

The resultant two pieces of diagrams represents the products of A(2,i)(x + (2j + l)z) 
and A(^2,2)/{i)ix — {2j + l)i) apart from normalization factors. Note that 2j = (modP) 

for j even and 2j = + (niodP) for j odd . Then from the expression (j^), one 
finds for the first piece, 

A(2,i)(x + (2j + 1)^) = 0(x + (2j - \)i)Tb,{x - ^^^i). 

The second piece, A(2,2)/(i)(a; — (2j + l)i), equals to (A(2,i)/(i)(a; + (2j + l)i))* by corollary 
Therefore, after taking account of proper normalization factors and using the property 
(Tb3(x))* = Tb.,{x*), one finds, 

T (^^ ^-^^ - /^2j+2,2i+i,2j)/(2,i)(x) jeven, 
2L + A [^2i+2,2i+i,2i)/(2,i)W ^odd. 

Thus the second relation is proved from (fH^ and □ 
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For L = 4k — l,4:k and 4^ + 1, we introduce 




< j < A; - 3 



j odd 



J even 



< J < A; - 3. 



(34) 



(35) 



while these do not show up for L = 4k + 2. 
By noticing 2 = 2U^: ~^ immedia 




(36) 



(37) 



These are desired relations. See the Y— system in the appendix 1X1 

The s/2 T— system is not itself closed. Here we find it connected to the s/3 like T— 
system through Tb^{x) = Td^{x). The latter system is, however, neither closed. In the 
next section, other functional relations are introduced, referred to as "magnon-like" , 
which turns out to "close" the total T— systems. 

6 The magnon-like t— system and the Dynkin hke 
diagrams 

The above transfer matrices are not enough to obtain sets of closed functional relations, 
which are necessary to derive TBA. We have to introduce further objects separately de- 
pending on four families. They are related to the magnon-like Y— system introduced in 
|25| 1^ . A magnon-like Y— system is associated to a diagram of the Dynkin type. Al- 
though they can be classified into the four families, the diagrams are all different depending 
on values of L. Consequently, the magnon-like Y system are all distinct for different L. 

Here the situation is slightly unified. The diagrams we have to deal with are "tails" of 
those in ISSHSD], independent of L, except for the dilute models. 

Before starting discussions on individual cases, we shall make remarks. A blank node 
in each Dynkin like diagram indicates its different character from the rests. To node a 
in the diagram, we associate t^m{x), where m takes an integer value from a finite set. 
Our magnon-like t— system is a set of functional relations among t'^ix)^ which should 
be described separately for four families. The most crucial observation below will be that 
t^m associated to a blank node is always expressed by (product of ) Tgg . 
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The dilute A^k_i model and the dilute ^4^+1 model share a same diagram but possess 
different t— system. The only "L" depending diagram for ^4^+2 {L = Ak + 2) remarkably 
coincides with the Dynkin diagram for Dk+2- 

Below we take x to be real. As usual, the symbol a h means that nodes a and h are 
adjacent on the diagram. 

6.1 The magnon-like t— system for A/^k-i 



6 




7 2 4 5 3 1 

Figure 8: Dynkin like diagram for A^k_i 
We define a magnon-like t— system associated to fig. |H1 

ttlA^)t^:!+A^) + \{t^^{x), a = l~6, m = l,2 (38) 

1 (a = l~6), = 0, (a = 3~6) 
t?{x). (39) 

This is not a set of closed functional relation. We, however, note its affinity to the Eq 
related T— system Indeed if ti\x) = 1 and t^\x) = tf\x) = 0, the above functional 
relation coincides with the T— system for level 3 RSOS model, corresponding to the coset 
iE^^^)2 X {E^^^)i/{E^^^)3. Our interpretation is as follows. The nontrivial t^i^ (x) makes the 
t^^^ (x) and tg^"* (x) non zero and connects the Eq level 3 T— system to others. Thus t^i^ (x) 
plays a role of "glue". This seems to be quite parallel to the s/3 like T system, which is 
connected to s/2 T system by nontrivial T^g. Here the situation is indeed the same; we 
will see below t^i\x) = Tb^- Tb^ thus plays a role in gluing the s/3 , the Eq and the s/2 
T— systems. 



2L + 4' 



tP(x± 



2L + 4' 
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We shall define "gauge invariant" functions for later discussions. 

a = 1 ~ 6 



t'i\x)t'i\x) 

t^2\x)t?\x) 

t?ix)tf\x) 



6.2 The magnon-like t— system for A 



4fc 



4 







1 



o 

7 

Figure 9: Dynkin like diagram for A/j^^ 
Associated to fig. El we propose a magnon-like t— system for this case, 

±)(x±^) = ttUx)t^:,l,{x) +l[t^^\x), a = l~6, m = l,2, 

= 1 (a = 1-6), tl^\x)=0 (a = 0,3,4), 
t^ix). 



L + 2' 



L + 2' 



Again, this is not a set of closed functional relation. In this case we note its similarity to 
the T— system for level 3 RSOS model, corresponding to the coset {A^^^)2 x {A^^^)i/ {A^^)-^ 
[HT)] . As in the preceding case, ti [x) will turn out to be equal to Ts^ix). Thus the 
interpretation is the same; t\ connects the T— system for level 3 A^ RSOS model , the 
sis and the s/2 T— systems. 
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Y functions for this case read, 



t^\x)tf{x)t?\x) 



tf\x)tf\x) 



(40) 
(41) 
(42) 



6.3 The magnon-like t— system for A 



4fc+l 



O- 

7 



Figure 10: Dynkin like diagram for A^k+i 

Although it shares the same diagram with the case Aik_i , a magnon-like t— system 
takes a slightly different form, 



tf^(x± 



2L + 4 
4i 



2L + 4 
^ ^ 2L + 4 



t^'\x± 



2L + 4 

i 



t?\x± 
ti'\x± 



2L + 4 

i 



2L + 4 

i 



2L + 4 



tt\x) 



XI 



3,4,5,6 



t?{x)+tf\x) 
t?\x)+t^P{x)tP{x) 



t!\x) 

ll4\x) 



a = 5, 6 



tf{x)tf\x)+t'i\x)t^^\x± 



tf\x)tf\x)+t)^>{x)t\'>{x± 



(5), 



.(2)/ 



2L + 4' 
2L + A' 



They are analogous to the T— system for level 3 RSOS model, corresponding to the 



coset {D^P)2 X {D\'')i/{D\'% [HE], which is clear by setting t\'' = t['' = t['' = 1 and 



(1) 



.(2) 
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4^'* ~ 4^'' ~ 0- ■'■^ ^^is case, will turn out to be a product of Tb.^, 4^'* ^-iid t^^^ are shown 
to be originated from the s/3 structure. Again, our interpretation is that non-trivial Tb-j 
glues the D4, the s/3 and the s/2 T— systems. 
Y functions are defined by, 

-i^fa) / \ 4 '*(-^) o r 

:= — , a = 3 ~ 6 



(1), , tf)(x)f)(x) 



4 ''(^)4 ''(^ 2L+4^) 



(2). , ._ 4'^(a^)4'^(^) 



4 ^(^)4 \^ ^ 2L+4^) 



fc-2V-^; •- (1) 



yR(x 



k-2\-^> ■- (2), . 



6.4 The magnon-like t— system for ^4^+2 



O- 



12 3 k-1 



k + 1 



Figure 11: Dynkin like diagram for ^4^+2 

Only in this case, the diagram depends on k. A magnon-like t— system for the fig. E] 
is defined as , 

tt\x±^) = 4:^i(^)4ti(^) + n^-(^)' « = l~fc + l m = l,2,3(43) 

4"^ = 1 (a = 1 ~ A; + 1), t'f\x) = (a = 2~A; + l) 
tf)(x±-^) = 4°'(x), tPix±j^^)=tf\x)tf\x). (44) 

In this case, the functional relations looks like the level 4 RSOS T— system for -D^+i, 
instead of level 3. 

The description of vj"''^ can de done in the most systematic way among 4 categories. 
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6.5 The solutions to magnon-like t— system 

The main finding in this report is that the solutions to the above magnon-hke t system 
are expressible in terms of fusion transfer matrices appearing in dilute Al models. We 
summarize the explicit solutions in appendix El The proof requires several steps and 
rather lengthy, part of which will be given in section ^1 Before going into such details, 
we would like draw attention to a general property of the solution and also its implication 
for the Y— system . Firstly, we need to explain Tk{x) appearing in the solutions for L = 
even. 



7 Kink transfer matrix 

The case studies on M4 5,M6,7 [23 imply the necessity of the introduction of objects 
which are not produced from boxes by (|T9|) . Generally, we find it possible to introduce 
the desired objects for the Al model , L even. 

The explicit procedure procedure to construct is as follows. Let L = 2£. and consider 
the model in regime 2, corresponding the perturbation 0i 2. Take 

as the top term (or the highest term). Next, one adds a lower term 

£ + 2. 0(x-|±f^)Q(x-fiiO 



U 0[X 



^ + 1 Q{x - j^i)Q{x + 

so that the singularity at x = Xj + -^i cancels due to the Bethe ansatz equation (jH)). 
Then we introduce another counter term so as to kill the singularity at x = Xj — ^i^i, 
and so on. 

One observes that this procedure closes by finite steps and a pole-free object, consists 
of 2^ + 1 terms, results. We refer to this as an eigenvalue of "kink transfer matrix" , 
Tk{x), although the corresponding transfer matrix is not yet found in the lattice model. 
For example, 

Tk\x) = W(p[x + -i)— rT: + w(p[x — 



5'Q{x- It) 5 ' Q(x - lt)Q{x + |z) 

^ 2., ^2 Q(x+jOQ(x-Jz) 8 Qix-jt^Qix-jt) 

+w d>(x + -i) — 1— — 1— + wd>(x z) — — ^ 

5 ^ g(x - lt)Q{x + It) 5' Q{x- lt)Q{x + t) 

, , Qix)Qix - It) ^ , s .^ Qjx + lt)Q{x + It) 

+ t)Q(x - z) + '^^'^ + 5'^ Q{x + lt)Q{x - z) 

2 Q{x-lt)Q{x + li) Q Q^a: + li)Q{x + lt) 

+W (pix — -tj— TTTTZ: (pix + -t^ 



5 ' g(x + h)Qix - U) 5 ' g(x + U)Q{x - U) 



5 

+w (pix — 



5 ^Q{x+U) 
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for the dilute Ag model. 

The numerical investigation on the largest eigenvalue indicates {Q{x*))* = Q{x), which 
leads to an important property (Tk{x*))* = Tk{x). 

For our argument, it is not important whether it is a "real " transfer matrix or not, but 
significant are its pole-free property as well as functional relations among the products of 
kink transfer matrices and T^ix). One finds, for instance, 

Tk{x ± ^^^0 = 2T4fc_i(x) 

in the case of the dilute A^k model. 

These relations are summarized as lemma [3 in appendix O 

As mentioned there, although the explicit proof has been done only for smaller values 
of k, we assume the following. 

Conjecture 1. Lemma^ holds for general k. 

The functional relations there play an important role in the proof of the magnon-like 
t— system for L = even. 



8 The reality property of magnon-like t— functions 

We remark the following significant property of our solutions to the T system. 

Lemma 2. Let T{x) be an either solution of the magnon-like t— system , the sl^ like T— 
system, or breather T system. Then, in "the largest eigenvalue sector" , 

T{x) = {T{x*)y (45) 

which we call the reality property. Namely it assures the reality of T{x) on the real axis. 

As is commented in the previous section, Tk{x) satisfies the reality property. One 
also verifies the property for most of tm^ (x) listed in the appendix El using the diagram- 
matic symmetry (PT|) . immediately. It can be shown for the remaining tm\x), which takes 
"asymmetric" shape , with utilizing the magnon-like "t-" system which can be established 
without use of the reality property. We take the dilute A^k case for example. Using the 
diagrammatic symmetry, tm\x) are shown to satisfy P5|l except for (a, m) = (0,1) and 
(a,m) = (4, 1). The one of the magnon-like t— system, 

t^l\x±^^^)-t^\x)=tf\x)tf\x) 

can be established independent of the property. Since all other entries other than tf'\x) 
satisfy , it follows from the above relation that tf'^ {x) also should satisfy . The 
case (a, m) = (4, 1) can be treated similarly. There are, however, two exceptions, t\ (x) 
and t^i\x) for the dilute ^4^+1 model. One still can devise other functional relations which 
explain the reality property of them. Since the proof is rather technical, we omit here. 
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9 The Y— system and the T— system 



Several Y— functions are already introduced for 4 categories by appropriate ratios of t 
or T functions. Thanks to the T— system or the t— system, one can easily check that 
the Y— system is satisfied for the case in which the Ihs are products of Yj or Yj'^\x) 
except for j = k — 2. One can verify the rest of the Y— system by utilizing the explicit 
expressions of t in appendix |B| 

Our main message in this report is summarized as follows. 

Proposition 1. The four categories of the Y— system proposed in 125] can he solved in 
terms of transfer matrices appeared in dilute Al model, if x is real. 

Due to the different forms of Y— system, the proof must be done separately for four 
categories. Below we shall concentrate on the case L = 4k . The other cases can be 
established with minor modifications. 

Proof. Proposition ^ for L = Ak 

Thanks to explicit relations between Y and t in ()40|) - ()42j) the "tail " part of the Y— 
system (j59l ) is shown to be a consequence of the magnon-like t system. Eq. (j58l ) also 
follows from the s/3 structure (jHEl ) except for j = /c — 2 . Consider the case a = 1 and 
j = k — 2. By replacing the Y functions by tm^ , one finds the Ihs, 



while the rhs reads 



^2 ''(a: ± j^^i) 
f)(x)t«(x±^z)(l + ^)- 

At this stage, the equality is not obvious. By the exphcit forms of yI^^{x) in ). we 
rewrite 

k odd 



T2(fc_2)(x±20 



1 _ I ^02(fc-2)(^) ' 

By substituting further the explicit forms of tm^ in appendix IB. 21 , one verifies the case 
a = 1 and j = — 2; for example both sides reduce to 



in case k even. The case a = 2 and j = k — 2 can be shown in the same manner. 
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Since eq. fl3Ul ) is already proved in , it remains to show eq. (fFFI ) . We must treat 
the cases k even and odd separately. For k = 2m + 1 , the product 

£=0 e=i 
is written in terms of fusion transfer matrices using (jH^ . 



n 



n 



£=l/=even 'C'2fc_2£'-2^'^ L+2 V ^ £'2fc-2£+2 '^*^ L+2 



n 



f^, -1- ii _|_ 4£'-l „-\7n ( I 4(^'-l)+l -x 



To(a;±ff^z)To(x±2'="5„ 



L+2 V-uv- ^ L+2 ^ 

where we have used T^^(x) = T^{x ± 2i) = To(x ± 4^0- Muhiplying this by the 
remaining term in the rhs, (1 + 1^3(3^)) = Tb^{x ± ^^)/To(x ± ^^^i), one reaches the 
expression identical to the Ihs. Note To^i^x) = Tb^{x). The case k even can be treated 
analogously. □ 

Once the solutions to magnon-like t system are proved to be given by those in the 
appendix El the solutions to Y system immediately follow. It is thus vital to verify the 
solutions. The proofs for the L = even cases and the odd cases are similar, respectively. 
We then present the proof for the case L = 4fc + 2 in the next section, where it can be 
done most systematically. As a representative of the odd case, we supplement the outline 
of the proof for L = Ak — 1 case in the appendix |El 



10 Proof of the magnon-like t— system for L = 4/c + 2 

We first prepare a few lemmas 
Lemma 3. 

2^77, -|- 5 

\m,m,l)/{m-l){x) = TD^_,{x)(f){x ± z) 

(2) 

This lemma holds not only for L = 4k + 2, but all L. The proof is direct by the a)^ 
property in (fTTl ). Although it is not used in this section, this lemma is frequently used in 
the proof of the magnon-like t— system. See appendix O for example. 
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Lemma 4. The following equality holds by mod P 



k 



2k + 2 
j 

2k + 2 



2k + 3, 
2fc + 3 + f, 



k 
k 



even 
odd 



2j, 

2j + ^ 



2 ' 



even 
odd. 



(46) 
(47) 



They can be checked trivially. 



Lemma 5. 



Lemma 6. 



tf\x)tt\x) 



2k+2 



[X 



k 
k 



even 
odd 



Do 



[x)), k : even 
[x)), k : odd 



(4J 



The proofs of the above two lemmas are given in the appendix iDl 
With these preparations, we shall prove the magnon-like t system equation by equation. 
Consider first the case m = 1 in 1)431 ). For a = 1, ■ ■ ■ , — 2, by paying attention to (fTTI 
) with j = 1, one immediately sees that the assertion is actually equivalent to (jHUl ). 



Thus only the cases a 
relation (jlflll ) with j 



k,k — 1 need proofs. Using t\ 
k, we find 



X) 



Tk{x)/2 and the functional 



t{'{x± 



1 



X] 



X] 



2k + 2 



T2k-2{x) + \{T2k+2{x) 



^2^1-2 (a;)) k '■ even 
T2k-2{x)) k : odd. 



The first and the second term in the rhs coincides with t^"" '''{x) and t^\x), respectively, 
which proves the case a = of (jlHl ). The a = k — 1 case is less trivial. Substituting 
a = k — 1 in and using (pll ) , one finds 



_^ 1 _ I T2k-2{x ± 2i) = r2fc(x)T2fc„4(x) + Td^^_^{x) keven 
' 2k + 2>-\T^,_,ix±2t) = T^,ix)T^,_,ix) + Tl(x) kodd. 



The second term in the rhs is identified with ^\x)- 
We note 



t['\x)t['+'\x) = {TK{x)/2f 




keven 
fcodd 



where (jlUll ) with j = + 1 is used in the second equality. Since T2k-A{x) = t^^ the 
first term in ()49|) can be written as t^i~'^\x)tf^\x)tf^~^^\x). Combining these observations. 
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we find that the rhs of coincides with {^) +4*^ {x)tf^^^^ {x) , which proves 

(jlSl) for a = A; - l,m = 1 . 

Next, consider the case m = 2 in PHl ). As in the case of m = 1 , the relation is 
equivalent to (jHTl ) for a = 1 ~ /c — 2. One only has to pay an attention, 

noix) = Toix ± ^^z) = t^^\x) 
for a = 1. In the case of a = A; — 1, the product of t^2 ^\^) reads. 



'n^^Jx ± 20 = T^,_2ix)TB,ix ± 2iT2^) + njx)n,,Jx) ^even 
2A; + 2 ' 1 Td^^^,,{x ± 2i) = T2k-2{x)TB,{x ± 2^i) + Td2^{x)Td,,,_^{x) kodd 



2k+2' 

where is again applied. The first term in the rhs is the product t[^~^\x)t^^~'^\x) , while 
the second is t'^\x)t^2^^\x)t'2~'^\x)- This comes from Lemma El as t^2\x)'t^2^^\x) = 
{t^\x)y. Thus the case a = k — l,m = 2 is proved. 

The case a = fc, m = 2 is the most non-trivial. We treat only the case k even for brevity. 



x). 



When taking the product (xiij^), we substitute two different expressions for t 

The first comes from the already proved relation, t^2\x) = t^i\x —t^i^^\x). The 

second is as listed in 

The first expression results, after applying (jlOlj) . 



t^2\x±i 



2k + 2' 



= ^(T2k{x)T2k-4{x) + T2k{x)T2k+4{x) 

1 1 1 1 

-'^2k-2{x - ^^^T^k+2{X + - T^k-2{X + ^^^)T2k+2{x - 

= -(T2k{x)T2k-4{x) + T2k{x)T2k+4{x) 

-T2k-2{X - 2t)T2k+2{x + 2l) - T2k-2{X + 2l)T2k+2{x - 2z)) . (50) 

From the second expression, t'^\x) = {T2i.j^2ix) — T2^_2(a;))/2, we find 



^ ^ 2k + 2' 4 



^{T2k{x)T2k-A{x) + T2k{x)T2k+i{x) + Td^^_^{x) + Td^,^^{x) 
-T2k-2{x - 2i)T2k+2{x + 2i)) - T2k-2{x + 2i)T2k+2{x - 2i)pl) 



where (jHIljl is used. Twice (jHT|) and subtract (jKHjl . one arrives at 
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The first term in the rhs coincides with ^3 (x) , while the second agrees [x] 
due to Lemma El Thus the case a = k,m = 2 is established. 

Finally consider the case m = 3. Except for the case a = 1 , the relation is satisfied 
trivially. When a = 1 , the Ihs reads 



(I) 1 k — 2 k 

t\'(x± — zi) = Tb^[x ± i)TB^{x ± 



2k + 2 ' 2A; + 2 ' 2A; + 2 ' 

k — 2 2 
Tbs{x ± ^^-^{}{To{x)To{x ± ^^r^O + TbAx)Tb.Ax)) 

k — 2 2 k — 2 

Ts^ix ± ^^-^{)To{x)To{x ± ^j^^^ i) + Td2{x)TbA^)TboSx ± ^ ^ 



where (I2SI) is used in the second quality and Td^^x) = Tb^{x) is used in the last. 

2 

2fc+2 



By noting tf\x) = To{x)Tq{x ± oFi^Oj find the resultant expression coincides 



with tf\x)tf\x)+t^^\x)t^^\x). 

We thereby verified the magnon-like t system for L = 4k + 2. 



11 Thermodynamic Bethe Ansatz in the scaHng hmit 
11.1 Analytic properties 

The Y— systems are obtained in |2S] by the simply Fourier transformation of the thermo- 
dynamic Bethe ansatz equation (TBA, for short), which is proposed in advance. In the 
present report, the y— system has been proved first, in the preceding sections. The (in- 
verse) transformation of the y— system to TBA seems thus to be direct. This is, however, 
not necessary true. One needs further analytic properties of of Y— functions as we shall 
see below. 

The F— system typically assumes the following form. (See appendix^for notations.) 

Fjr) {x±za) = l[ Sf ) {x ± ^^f^ ) J] 4'^ ± ) (52) 

b,j c,k 

We try to take the logarithm of both sides and take the Fourier transformation. For the 
necessary condition that the last transformation can be done simply, one must ensure 
that the Ihs (rhs) is analytic and nonzero in the strip Qx G [—«,«] {Qx G )• 
Otherwise, one must take account of contributions from zeros in the strips. 

As is already shown, complex combinations of "transfer matrices" solve the both sides. 
The explicit expressions for eigenvalues of transfer matrices thus enable the investigations 
of the analytic properties of the Y— systems, once parameters, N, q, (3 are fixed. We have 
carried out numerical investigation for ranges of parameters, L = 5,7, 8, 9, 10, 4 < < 
10, < g < 0.5, f3 = 0.04, 0, 05, ■ ■ ■ , 0.1. Although our main interest is focus on the limit 
A" — *• 00, our data suggests weak dependency on A^. For illustration, appendix IHl presents 
some numerical data and graph showing locations of zeros for various transfer matrices 
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for the dilute Ag model in regime 2 . The analytic property of Y can be easily read off 
from these zeros. 

To state our conclusion, we remark that the Y— system is invariant , for even A^, if Y 
is replaced by Y, defined by 

I . ^■7^^#^^ for u<0 



Ybj^{x)k{x ± i{l — '»0 2(L+^2) ) fo^ M>0 



and all other cases, Y = Y . The parameter v! stands for ^^i^m- This is due to the 
definition of 



K\X) — \% fL+2~) 

^2(^i(izi)7rx,r')- 

which satisfies k{x ± i ^[^2) ) ~ '^^^ elliptic nome g' = exp(— r'),r' = ^^z^t is 
introduced so as to respect the periodicity of the Y function on the real direction of x. 
We refer the resultant functional relation to Y- system. 
Our conjecture is as follows. 

Conjecture 2. The both sides of the Y -system are analytic and nonzero for appropriate 
strips, for any L, N,q, p. 

11.2 The coupled Integral equations 

Once the conjecture |2l is taken for granted, the transformation of the F-system into a set 
of integral equations is straightforward. When all Y functions is replaced by Y, only the 
integral equation for log Yb^ (x) possesses a drive term. That is, for m > 0, with appropriate 
kernel functions KBi,a,m-, it takes of the form, 

L — 2 1"^" 

^ogYs.ix) = -\ogK{x±i{l-u')—j-—-)+^ KB,,a,m{x - x')\ogE^^\x')dx' 

{a,m) 

+ J^B,,aA^-x)\ogCt\x')dx' (53) 



{a,m) 

where r" = ^. For u < 0, the driving term should be logfi;(x ± i{l + u') 2(1+2) '^ ■ 
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The free-energy is evaluated from 

logTsiix) = logrN{x) + J s{x — x')logeNix')dx' 

+ J s{x~ x')\og{l + YB,ix'))dx' (54) 

/ h ( nr* - 

rN{x) 

ENix) 



h{x±C^,+u)t)h{x±{^,+u)t^~ ^ 



h{x±{^,-u)mx±C^,-u)i] - - 
Mx±(^ + «)z)Mx±(|f±| + «)z) 
1 



2 cosh 2^ A;' 



Here we denote the Fourier transformation of F{x) by F{k), 

F{k) = [ F{x)e-'^''dx 



F{x) = l-Y^F{K)e'^-\ K = n5 

n 

and S = vr/r". The free energy is determined by the value of Tbi{x) at the origin, (3f = 
— logTB^(O), thus we are only interested in x ~ in (jSH). 

All the above results are valid for any A^. Now the Trotter limit can be performed 
analytically. The drive term in ()53|) becomes, 

L — 2 

- log k{x ± z(l - ^') 2(^i_^2) ^ ^ 4/37rs(x) 

while 

/ 1—2 21 

logr^(O) A/3((log^i(-^A))'+(log^i(-^A)y 

ioi u = (3/N . The crossing parameter A is set to be f|^^- 
11.3 Thermodynamic Bethe Ansatz equations 

We consider the scaling limit g ^ 0,/5 — *• oo ^H]- The precise tuning of these two pa- 
rameters will be specified below. In order to take this limit, we conveniently rewrite the 
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integrals in the preceding subsection in such a form that "kernel function times small 
functions" . The numerical calculation indicates, in the vicinity of the origin, 




if /? >> 1. In the connection with Ml^l+i + $1^2 theory, we are interested in the regime 
2 (m > 0) case. Thus we first rewrite the Fsystem such that the relations involves only 
(1 + y). Now the important contribution for the integral comes from "Fermi" surface 
x' ~ r". Thus we are interested in the scaling function, \img_,Q yBi{x + T") etc and integral 
equations among them, TBA equation. To derive them, we take the logarithm of both 
sides of the transformed Y system above, take also the Fourier transform and find 



(55) 



where Lb^ = log(l + and similarly for others. Note that M and Kq are asymmetric 

matrices. See appendix El for the example for L = 5,7. By multiplying from the 
left, @the kernel matrix for TBA, M~^Kq turns out to be symmetric, remarkably. This 
property is crucial in applying the dilogarithm technique to evaluate the central charge. 
Its elements agree with the expression described in |25j in terms of S matrices, under 
identification x = 36 /tt in the limit q —>■ 0. Finally, we have checked that nonzero drive 
terms exist only for Y functions corresponding to Breathers and Kinks. For example, we 
consider L = 5. The explicit forms of the drive terms in the Fourier space read. 





/loj 






(l\ 






M 


loj 















V 


'■■ J 




v) 




V '■■ J 



d 



Bi 



SvT/S cosh jifc 
(2cosh^A;- 1)D(A;) 



d 



B3 



IB5 



14 

47r/5(2cosh^A; 



1) (2 cosh ^ A; 



d 



B2 



''Ki 



d 



K2 



D{k) 

IGnP cosh jgk cosh j^k 
{2cosh^k-l)D{k) 

SnP cosh jjk 
{2cosh^k-l)D{k) 

AttP 

(2cosh^A;- 

SttP cosh Y^k 
(2 cosh Y^k 



i)D{ky 



where we denote by dsi for the drive term associated to logY^j and so on. A common 
denominator D{k) denotes 



Dik) = 2 cosh —k + 2 cosh —k 
^ ^ 14 14 



6 4 
2 cosh —k — 2 cosh —k + 1. 
14 14 



33 



When taking the inverse Fourier transformation, the nearest zero to the real axis of the 
denominator is relevant in the "scaling" limit, q —* 0. We have checked for various L that 
the zero of the common denominator D{k) brings the nearest zero, and always, k = ±^|. 

Let us be more explicit for L = 5. As announced above, we are interested in the 
integral equation for y{x) = limg^o Y{x + r"). Thus the associated drive term reads, 

^scaiing^^^) := lim d{x + r") = lim — y / 6^''''^^+''-''+^"^''' d{k')dk' 

where the Poisson's summation formula is applied. As g — > 0, 6, r" tends to be infinity, 
the only m = 0, m = — 1 terms contribute, 

^scaling _ / / q-'^^" -''U{k')dk') . 

q^0\2TTj_^ 2tcJ_^ J 

We enclose the complex contour of the integration in the upper half plane for the first 
term and in the lower half plane for the second. For r" oo, the nearest poles k = ±z7r/3 
bring the dominant contributions. We denote 

1 r 



D{k) {k-i^)i' 



3' 

then 

jscaling/ n S'^'" -±t i ^ 

"k- (X) = e 7 cosh — X 

^ ' 2cos^-l 3 

for example. Let R be radius of compactification. Then we take q ^ 0, (3 ^ oo while 
mfc finite; mkR = 2^^^^^- This is the precise meaning of our scaling limit. We write 
|x = 6. Immediately, one verifies that all other drive terms also take the form mi? cosh ^ 
and their mass ratio agree with those in |25j . 

11 . 11 3 

ruB^ = 2mK cos —TT = 4mif cos — ttcos — tt 

1 4 1 

niBs = 4:mK cos — vr cos — vr = ^uik cos —n 

4 

= 2mK cos 



We then write ya{6) = exp(e(0)), and verify that the TBA is recovered in the scaling 
limit. 

2{L+1) 

We comment that for general L, the straightforward calculation shows mx oc g3(L+2) _ 
This is consistent with scaling relation, ^ = ^ for m oc q^. Indeed comparison of the 

result leads to proper scaling dimension, A 



il,2 

L-2 



4{L+1) • 

It remains to check the expression for the scaling free energy. 
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The scahng hmit of free energy comes from the contribution from the last term in 
) at X = 0. 

/oo /"OO /"OO 

s{x')\og{l + YB,{x'))dx' = / s{x')\ogYB,{x')dx' + / s{x') log LB.dx'. 
■oo J —oo J —OO 

We substitute the inverse transformed (fS^ into the first term in the rhs above to find, 

/oo roo poo 

s(x') log(l + YB,ix'))dx' = / s{x')dB^ix')dx' + / }Ce{x') log L(,dx'. 
-oo J —oo £ J —oo 

where JCi reads in the Fourier space, 
and 

A,,- := M,j + {Ko)i,j. 

We adopt indices i/Bi = Ui^Ubz = 2/2, and so on. 

The scahng free energy /s, the contribution of excitations near the Fermi surface, is 
thus identified with 

1 r°° 1 

/, = -hm-V / /C,(x' + r")log(l + -(x'))rfa;'. 

3-^0 p ^ J-oo yt 

To perform the hmit, we use a remarkable relation, 

K,i = de 

which we checked case by case. Then one repeats the same argument in deriving the drive 
terms in TBA equation, and reaches 

1 /*oo -1 -1 roo 

fs = y" / m^cosh^log(l + - )c/^,= y/ mfCosh^log(l + e-^(^))(i^ 

under tuning /3 = ^J^R. Therefore the field theoretical expression of the free energy is 
also reproduced. 

12 Summary and discussions 

In this report, we establish relevant functional relations of the dilute models for ar- 
bitrary L, to analyze finite temperature property we demonstrate explicitly that TBA 
equations for Ml,l+i + 0i,2 , conjectured by Dorey, Pocklington and Tateo [23], is realized 
in the scaling limit. The crucial idea is to introduce fusion transfer matrices associated to 
skew Young tableaux and to investigate the functional relations among them. 
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There are still many open problems. The explicit identification of string solutions would 
be definitely one of the most important. The complete study on this will shed some light 
on the way how to proceed for TBA in the case of perturbed nonunitary minimal models. 
We mention the first step in this direction in [nZj. The minimal unitary theory perturbed 
by (j)2,i may be treated within the same framework. One only has to deal with the regime 
1 instead of the regime 2. We already have some results about the "exceptional" cases 
in terminology of [23] , however, they are omitted in this report for brevity. We hope to 
complete the program on the 02,1 theory for arbitrary L in a subsequent report. 

Before closing this paper, we comment on a possible link to Fermionic formulae 
of the Virasoro characters. In EHl lEHl IMl IBS], remarkable results are reported on 
new representations of the Virasoro characters which originate form fermionic parti- 
cle bases of the Bethe ansatz equations. See also combinatorial arguments in view of 
TBA |63| IMj and mathematical generalizations based on the combinatorial objects, the 
rigged configurations [HKj- . The Fermionic characters obtained in [SHI seem to 

be connected to (^1,3 perturbation. In [SIl EII, they also explicitly showed the different 
realization of Fermionic formulae which are related to 02,i or 0i,5 perturbations. 

Since the TBA in offers a description of the associated Hilbert space, it might 
be natural to expect that a new fermionic representation is possible based on the 0i,2 
perturbation. There exists already positive results based on the simplest case, the dilute 
A3 model [72]. It may be natural to expect the generalizations. 

Indeed, we check up to 0(g^°°) that the following formula is valid, 

2 

where 



9 = 





6, 


8, 


4, 


2, 


4, 


0, 


o\ 


6, 


12, 


16, 


8, 


4, 


8, 


0, 





8, 


16, 


24, 


12, 


6, 


12, 


0, 





4, 


8, 


12, 


8, 


4, 


6, 


0, 





2, 


4, 


6, 


4, 


3, 


3, 


1 

2' 





4, 


8, 


12, 


6, 


3, 


r. 


1 

2' 





0, 


0, 


0, 


0, 


1 

2' 


1 

2' 


1, 


1 

2 


^0, 


0, 


0, 


0, 


0, 


0, 


1 

2' 


1 / 



The summation is taken under the condition, < "5+'^6+"8 ^ Z>Q,ns < ^ G Z^q. Note 
that g has its origin in the kernel functions of the TBA equation. 

A similar formula for Xii is checked up to 0(g^°°), which may imply the existence of 
Fermionic formulae for L odd case in general . We are not yet able to find such a simple 
expression for the L even case. We hope to clarify this in a future publication ^. 
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^We have learnt that S.O. Warnaar has independently found polynomial version of these formulas 
corresponding to both $1,2 and <I>2,i perturbations. 
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A Y-system 



As obtained in j25], Y system for the minimal unitary theories perturbed by $1^2 falls 
into 4 categories. Below we shall summarize their result with slight changes in notations. 
Instead of their rapidity variable 6, we use x as in the text. They are simply related by 



X 



^. By C we mean C = for a Y function. In the section an inverse function 

L = 1 + y is also introduced but it will not be utilized here. We understand indices 
associated with C specify these for Y functions, e.g.. 



etc. Similarly Em denotes 



1 



XI 



^^:^\x) = i+Yi^\x) 



A.l Case, (L = 4A; - > 2): 

The Y-system is 



^ ^ 2L + 4 



X] 



k-2 



1=0 

k-2 



2L + 4 



n ^ 



-I s 



2L + 4 



2L + 4 



2L + 4 



'j,fc-2<'al 



X 



Cn,{x)CfUx)4Ux ± ^^)C^U- ± ± -^^) 



-I s 



2L + A' 2L + 4 

(with j = 0, . . . , k — 2] a = 1,2 and a = 3 — a) 



2L + 4 



1 



2L + 4 



nK-i(^)] (7=1,... 6) 
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where the last product is taken over nodes on the Dynkin diagram, which is obtained 
from fig.^lby deleting the node 7. /["^^-i] stands for the incidence matrix for A^-i, where 
nodes are indexed by to A; — 2. 



A.2 Case, (L = 4A;,/c > 2): 
The Y-system is 



L-2 , 
4 



2L + 4 



2L + 4 



£=0 



2L + 4 



1=0 



(1) 



2L + 4 



(with j = 0, . . . , — 2; a = 1, 2 and a = 3 — a) (58) 
''''^'''nK-i(^)] (7 = 0,... 4) (59) 



(2) 



'5i,fc-2'5cil 



(56) 
(57) 



where the last product is taken over nodes on the Dynkin diagram, which is obtained 
from fig. El by deleting the node 7. 
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A.3 Case, (L = 4A;+ > 2): 
The Y-system is 



4 

2L + 4 



2L + 4 
3 

2L + A 
1 

2L + 4' 
1 

2L + 4' 
1 

2L + 4' 
1 



^^3 



k-2 



5..M5£.WnHf'(x± '"-''^-^ l 



£=0 



2L + 4 



fc-2 



^=0 



X 



X 



2L + 4 

(with j = 0, . . . , k — 2; a = 1,2 and a = 3 — a) 



2L + 4 



2L + 4 



2L + 4 



'^fc-i("'')'^fc-i(-^) 



2L + 4 

A.4 Case, (L = 4A; + 2,A; > 2): 



L-2 , 
L-2 , 

^^^(^^2m'^ 



fc-i 



^=0 



2L + 4 



(with j = 0, . . . , /c; a = 1,2 and a = 3 — a) 

where the last product is taken over nodes on the Dk+i Dynkin diagram, which is 
obtained from fig. HUby deleting the node 0. 
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B explicit solutions to magnon-like t— system 



Below, sets of explicit solutions to magnon-like t— system are listed. There are, however, 
many equivalent expressions as remarked in the text. We write just representatives among 
them. When other expressions are convenient, some remarked will be made then. 



B.l solutions for A 



tf\x 



tf\x 



t^\x 



X 



2k\ 



k even 
k odd 

k even 
k odd 



A(6fe-2,2fc)/(l)(a^ 



2L+4^' 
2L+4 V' 



k even 
k odd 



2 2L+4 



)^*-(4fc,4A:,2A;)/{4fc-l,l) 



ii\ix — 2ki + i 



2L+4/' 



A 



(4fc,4fc,2fc)/(4fc-l,l) 



'x — 2ki + i 



2L + 4:' 



k even 
k odd 



A(10fc-4,6fc-2,6fc-3)/(4A:-l,4A:-2) (3;) , 
A(10fc-4,6fc-2,6fc-3)/(4fc-l,4fe-2)(^)' 



[x-{L + 2)i 



2i 



-A 



k even 
k odd 

L + 4: 



2L+4 



) 



{L,l)( 



2L + 4 



For a guide to eye, corresponding diagrams are illustrated in fig. [121 for k = 2. 



n 



(60) 
(61) 
(62) 

(63) 

(64) 

(65) 
(66) 











Figure 12: Young diagrams associated to t^i\x) for k = 2 
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For , 





[T2k-2{ 


t?\x) = < 




t?\x) = < 


TD2k-2 
/D2k-2 


tf\x) = < 


rp\J 
/D2k-2 


t^Pix) = t 


f\x± 


tf{x) = < 


'^D2k-2 
P^D2k-2 



2L+4 
X± ^ 



2L+4 V) 
X ± TTT^i), 



[x± 



2L+4' 
3 

2L+4' 



k even 
/c odd 

k even 
A; odd 



2L+4 
2 



2L+4 



^)^2^-2(^)' 



1 



2L + 4 



2L+4' 

i) 



k even 
A; odd 



k even 
/c odd 

k even 
A; odd 



and 



tf{x) 




k even 
A; odd 

k even 
A; odd. 



B.2 solutions for A 



tf\x 



t?{x 



t?\x 
t?ix 



t?{x 



tf\x 



ns ror /±L=4k 

'tf{x - ^)T2k{x + ^) - 4'^(^ + ^)T^k- 
t?{x - ^)TV (x + ^) - t?{x + ^)T,,_ 



i(x - j^),keven 
i'(^+z^)7^2fc-2(x-;^),A;odd 



|7^2".(^), 

I TD2k{x), 



k even 
A; odd 

k even 
/c odd 



0(x - ^i) 

= Tb,{x) 



L + 2 ' 2 
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Among t^\x), non-trivial are the cases a — 3, A. 



T2k-2{x), A; odd 

TD^k-2{^)^ /seven 
n,,Jx), A; odd 



and others are easily obtained from them due to the t— system, 

tf\x) = t^^\x)t^^\x) 
t?ix) = t^^\x±^) 

t?\x) = t^^\x± 



L + 2' 

i 



L + 2' 



Finally, 



{Tl_^{x)) T2k-4{x), keven 
{T2u-2{x)YtI_,{x), A; odd 



[X) 



k even 



4 odd. 



B.3 solutions for Al= 



■■ik+\ 



^2fc-2(^)' 
T2k-2{X), 



T, 



Do 



A(6fc+2,2fc)/(l)(a: 

-^(6it+2,2fe)/(l) 
1 



k even 
k odd 

A; even 
A; odd 

2L+4^)' 



2L+4' 



-A 



A; even 
A; odd 



[x ■ 



(4A; + 3)z) 



-i,i)(^+9) 



A 



(10fc+2,6fc+2,6fc+l) /(4fc+l,4fc) 



X 



■^(10fc+2,6fe+2,6fc+l)/(4fc+l,4fc) {x) , 



A; even 
A; odd 

3L 
2L + 4 



rB3(x± 



3i 
2L + 4 



) 



42 



O 



B 



Q 







Tl 











Figure 13: Young diagrams associated to t^i\x) for = 2 



The corresponding diagrams are illustrated in fig. El for k = 2. Notice their similarity 
to those in fig. [T21 



For t 



(a) 
2 



.(4), 



X) 



t^^Hx) 



t^^Hx) 



and the rests are found to be 



T2k{x)T2(k-2){x) 



D. 



njx)Tr 



2(fc-2) 



[X 



n,^{x)T2u{x± 2 



2L+4 

2 „ 



fx ± 7rF^0^2fc(a;) 



2L+4 



2L+4 



2k\ 



k even 
A; odd 

k even 
A; odd 

), k even 
), A; odd 

k even 
/c odd 



2L + 4 



njx)T2k{x), 



{x)T, 



2k 



X) 



k even 
k odd 



tf)(x)T53(x) 



(84) 
(85) 
(86) 
(87) 



^9) 



(90) 
(91) 
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B.4 solutions for AL^4k+2 

In this case, the eigenvalues of t associated to the tails, k and A; + 1 , are identical. 

t['\x) = tf^'\x) = ^^ (92) 

ja) ^ iT2a{x), a even < a < A; - 1 

\rv(x), a odd ^ ^ 

[TD2a{x), «odd 

(96) 

tf(a;) = tt'\x) = Tb,{x) (97) 
4«)(a;) = ^^3(^±^0, l<a<k-l (98) 

tj)(a:) = Ts,(a;)Ts3(a;±^^z). (99) 

C Functional relations for the kink transfer matrix 

We verify the following lemma for relatively small values of k 

Lemma 7. The following functional relations among kink transfer matrices and are 
valid for the dilute A^^ and A4k+2 models (k — 1,2,3 for the former and A; = ~ 3 for 
the latter) in regime 2. (corresponding to M^^^+i + 4)1,2) 
A^k ■ 

"^-^^^ ^ ilr^^^ " 2T4fe_i(x) 

T,U±^±_A±1,) = f2(T,,_,(x)+rv_,.(x)), j: odd ^.^1... ^ 

^ 4A; + 2 \2(rv._,(x)+r4.-2,(x)), J : even' ^ ' 

2(fc + j) + l _ j2(T2,_i(x)+T4V2,-i(^)), j: odd ._ 

Ak + 2 \2m_,{x)+T,,.2j-i{x)), j: even ^ ''■"''^ 
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A, 



4fc+2 



Tk{x±^^i)=2TI^,{x) (100) 



Tk{x± 



2k + 2 

k + l- 3 .N ^ \2{T2j^2{,x) + T4fc_2j+2(a;)), j 
Tf , fc + 1 + j .X _ /2(T2j„i(x) + T4fc_2j+i(x)), j 



odd 
even 

odd 
even 



j = l,---,A; + l (101) 
j = l,---,A; + l (102) 



We also note another functional relation for Tk{x ± i) of the case, indicating the 

03 property 




'^^ ^'■^^x)TKix) + 2T2k-iix)), k : even. 



There are also relations obtained from the above by using the duality and the periodicity, 
omitted for brevity. 

Proof. The proof is done by comparing the exphcit dressed vacuum forms in the both sides. 
In most cases, this procedure , although simple enough in disguise , is not straightforward 
. For example, for the dilute Aiq model, the product of two Tk with identical spectral 
parameters, decomposes as 

T^{x) = 3T,ix)+T,,ix). 

By the duality, Ti^{x) = T^i^x), one concludes p01|) . with k = 2 and j = 3. Similarly, 
consider 

T|(x ±'-)= 2T2{x + ^i) + Te{x + ^z) + T,s{x + ^z). 



The rhs is identical to 2(T2(x + yi) + Tq{x + yz)) due to the duality, which proves (HOH), 
with k = 2 and j = 2. The rests are also shown by using this kind of trick. □ 

Corollary 2. For Aik+2 models, the kink transfer matrix is represented as 

Tk{x) = < ^ ^'"'^^ " (104) 

l^^(^2fc+i(a;) - T2k-i{x)), k : even. 

This is shown by comparing ()102j) with j = k + 1 and p()3|l . 

For the dilute Al models with odd L, we can also construct Tk{x) following the above 
manner. The resulting pole-free object consists of 2L + 2 terms. It, however, proves to be 
identically zero under some assumption on analyticity. At the present, we do not know 
what is the implication of this. 
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D Proofs of Lemma El and M 

Proof. : Lemma El 

We use one of the magnon-like t system 

which has been shown in the main text. Instead of using the solution for t^2 \^) )• 
we solve the above for t^\x) and take its square, 

+ {tr'>{x))'' 

where the solutions for t^^'' is substituted. By further substituting the solutions for t^^~^^ 
and making use of the functional relations for ( the cases j = k and j = A; + 1 in 
eq. fjlOll ))■ one arrives 



(tr(x)) 



T2k{x ± jj^i) - T2^fc+i)(a;)T2^(fe_i)(x), k : even 
T^f^{x ± 2^z) - T2(fc+i)(a;)T2(fc_i)(x), /c : odd. 



By noting the mod P relation ()47|1 with j = 1 and using the s/3 type relation (jSDI), 
one concludes {t^\x)Y = T^^^{x) for /c even and {t^\x)Y = To^ki^) for k odd, which 
proves lemma |31 □ 

Proof. : Lemma ini 

Thanks to the (jr" property p()4j] . t^i\x) has an expression, 

^{k),, ^ Tk{x) ^ f 2^(T2fc+i(a;) - T2k-i{x)), k : even 
2 \24^^{T^k+ii^)-T,\_,ix)), A:: odd. 

We solve for t^^\x) = Tb-^{x) and substitute the result into the product , 

2i<t>^(Ti{x ± ^i)T2k+i{x) -Ti{x± ^i)T2k-l{x) 

^^H^^)-^2A:+i(a^) + Tq{x ± i^^i)T2k-i{x)) k : even 



tf\x)tt\x) 



2((^v^{a;))2 -'- 2A:+2^)-^2^+l('^) "^li.^ 2k+2^)'^2k-li^) 

-Ux ± z|±|^)TV^^(x) + To(x ± ^^^^)T^^_^{x)) k : odd. 

(105) 
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We will show that the above expression coincide with the rhs of ()48|). For simplicity, only 
the case k even is considered below. Consider the second term in the rhs of ()105p first. 



Noting the expression Ti{x ± gFi^' 
k 



Ti(x± 



2k + 2 



i)T2k-i(x) =To(x± 



Ti{x ± 2ki) due to (jTfj) with m = k, we have 

-i>T2k+l{X) + A(^2k-l,2k-l,l)/{2k-2)[X) 



2k + 2 

, f2k{x - i) f2k{x + i) 
+ ~f 7Zr^(^k,2k)/i2k--i){x) + -— A(2fc,i)(x). (106) 

j2k-l[X) j2k-l[X) 



The graphical representation in fig. will be intuitively helpful. The second term in the 
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Figure 14: The graphical representation of eq. ()106p for A; = 3 . A figure a in a box specifies 
its spectral parameter x + ia. 



rhs of 



is identical to {(j)^ {x)YT]:,^^_^{x) due to LemmaElwith m = 2k — 1. The third 



and the fourth term in the rhs of p06l ) can de further decomposed with the aid of the 
quantum Jacobi-Trudi formula and mod P relations in Lemma |3] 

ro(x - {2k - l)i)T2k{x - i)Ti{x + 2ki) - To{x ± {2k - l)t)T2k+i{x) 

k + 2 k + 2 

Tq{x + , i)T2k{x - i)Ti{x - {2k + 3)i) - Tq{x ± ^, , J )T2k+i{x) 



2k + 2 



2k + 2 



for the third, and 

To{x + {2k + l)i)T2k{x + i)Ti{x - 2ki) - Tq{x ± {2k - l)i)T2k+i{x) 

k + 2 k + 2 

To{x - ^, , i)T2k{x + i)Ti{x + {2k + 3)i) - To(x ± ^, , J )T2k+i{x) 



2k + 2 



2k + 2 



for the fourth in ()106|1 . Next we consider the first term in the rhs of fllOBj) . To obtain a 

(2) 

non trivial decomposition, we remark the following expression for T2k+i{x) due to the 
property, 

T2k+l{x) = ^(j,^ 4fc+5-')^(2fc+1.2fc+l)(a^)- 

We also use the expression Ti{x ± ■^tpl'^) ~ Ti{x±. {2k + ?>)i) due to PHji . By these tricks, 
three rectangles now set in the right position so that non trivial decomposition occurs. 

^^^^^ 2k + 2 ^^^^^^^^^^ ^ </'(^±^^^O^i32fc+2(a;)+^o(x±^^i)T2fc_i(2;) 

+ncA^2T+i,2.)(^ + 2^) + <A^-+i,2.+i)/(i)(^ - 2^. (107) 
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The scalar function of the first term in the rhs can be also written as (j){x ± ^^^i) = 
The graphical representation in fig .Elwill be again helpful. The normalization factor 
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Figure 15: The graphical representation of eq. ()107l ) for A; = 2. A figure a in a box specifies 
its spectral parameter x + ia. 

ric in the third and the fourth term in p()7l ) reads, 

03(a; _ ,(2A; + 1)) 



rir 



f2k+i{x) Uf=i' Mx - 2(fc + 1)^ + 2jt) 



The upper index "non" in the third and the fourth term indicates the lack of their nor- 
malization factors. See for 02- The third term in ()1()7|1 is further transformed with 
aid of the relation, 

^(2fc+l,2fc)(^ + 2«) = ^'{2k,2k)ix + i)Tl{x + {2k + 3)2) — A"2fe_i,2A:-l,2A:-l)/(2fc-2)(^ + 3«) 

2k 

= f2k{x + i)Y[Mx- 2ki + 2ji)T^{x + {2k + ?,)i)T2k{x + i) 



2k-l 



-f2k-i{x)Mx + {2k + II Mx - 2ki + 2ji)T2k-i{x) 

i=i 

where we have used tableaux rule in the second equality. After taking account of nor- 
malization factors and mod P relations in Lemma |3] , the third term in ()107|) is given 

by 

k -\- 2 k + 2 

Tq{x - ^, , i)T2k{x + i)Ti{x + {2k + 3)i) - Tq{x ± ——-i)T2k-i{x). 



2k + 2 

Similarly the fourth term reads 
k + 2 



2k + 2 



To{x + 



i)T2k{x - i)Ti{x - {2k + 3)z) - To{x ± 



k + 2 



2k + 2 ' ' ' ' ' ' ' 2k + 2 

Therefore the first two terms in the rhs of ()105l ) add up to 

k k 
Ti{x ± ^, , i)T2k+i{x) -Ti{x± ^, , J )T2k~i{x) 



i)T2k-i{x). 



2k + 2 



2k + 2 



k + 2 



^x)) VJ^D^u+A^} -Td^^_^{x)) + To{x ± '] ' ^ i){T2k+i{x) - T2k-i{x)) 



2k + 2 
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Obviously the last two terms in the above cancel the third and the fourth term in (jlOSj) . 
Thus we arrive 

1 

2 

The case for odd k is done in parallel. 

□ 



E Proof of the magnon-like t system for L = 4/c + 1 



We summarize the proofs for fort he magnon-like t system, L = 4k + 1. Most of relations 

2L+4' 

are omitted except for a = 3, 4 . 
Note the useful mod P relations 



in (jHHj) . of which Ihs takes the form t^2\^^ 27+4^)' shown without difficulty, thus they 



4 2 P 

= -^ = 2-- (110) 
2L + 4 Ak + 3 2 ^ ^ 

which will be used frequently below. 

E.l the proof for t^^\x ± = tf\x) + t'i\x) 

Thanks to ()11()|1 . the Ihs equals to T2(fc_i)(x ± 2i) or T^^f,_^.^{x ± 2i), depending on k even 
or odd. Then the sis relation (j30|) proves the assertion. 

E.2 the proof for tf \x ± = t''^\x)tf\x) + tf\x) 

The proof is similar. One only has to apply (jSJ)- 

E.3 the proof for tf\x ± = t^^\x)tf\x) + tf\x) 

We consider the case k odd. Using the realness of t\ (x), the mod P relation ()1()8|1 and 
the diagrammatic symmetry, we rewrite the Ihs as 

A(6fc+2,2fc)/(i)(a; + {4:k + l)i)A(6fc+2,6fc+i)/(4A:+2)(a; - (4A; + 

while , thank to the duality, the second term in the rhs is given by 

tf\x) = Tn,,ix)nk+iix ± (Ak + 3)i). 
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Then the quantum Jacobi Trudi formula leads to 

tf\x±—^t)-tf\x) = 
IV 2L + 4 ^ ^ ^ ' 



T2{k-i){.x) \T2k+2{.x)TQk+i{x ± (4/c + 3)i) + T2fc_2(x)T6fc+3(x ± (4A; + 

- T2k{x - 2z)T6fc+i(x + {Ak + 3)t)nk+3{x - {Ak + 

- T2k{x + 2t)nk+i{x - {Ak + 3)t)nk+3{x + {Ak + l)z)) 

= T2(k~l)ix)A(^iok+2,6k+2,6k+l)/{4:k+l,4k){x). 

One can check the rhs of the above equation equals to t^i\x)t^i\x) . 

E.4 the proof for tf\x±^i) = t^^\x ± j^)tf\x) + tf\x)tf\x) 
Again we treat the case k odd. Since t'^\x) = tf\x), the equality to prove is, 

± = tf> (x) (tf ) (x) + t? {x ± ^^)) . (Ill) 

Substitute (IHUj) and use (plj) for t^^^ (x) , the Ihs is given by 



4'^(a;± 



1 



2^^40 = tf\x){T,,{x + {Ak + l)t)T,,ix + 2t))l^^__^^ 
= t'i\x)(T2k{x + {Ak + l)t)nk+i{x + {8k + 6)t 



2L+4 



where the duality is used in the second equality. By the (jTHl), it is readily shown that the 
rhs of the bracket above coincide with 

4'^(^) {^^i6k+2,2k)/il)i^ - JJ^) + T2k-2{X ± ^^^)) • 

It is immediate to check that the above equation coincides with the rhs of (jlllj) 
E.5 the proof for tf\x ± = tf\x) + tf\x) 

The first term in the rhs, tf\x) is identified with A(iofc+2,6fc+2,6fe+i)/(4fc+i,4fc)(a;) for k odd or 
^aoA:+2,6fc+2,6fc+i)/{4fc+i,4fc)(3^) k evcu. Wc instead consider 7^iofc+2,6fc+2,6fc+i)/{4fc+i,4fe)(a;). 
See (EUjl for their difference. Due to the tableaux rule, it clearly decomposes into three 
parts. The "height-three" part reduces to the shifted product of 03 (x) in (|T^ . The re- 
maining two parts are identified with the diagram for tf^ and its 180° rotation, referring 
to (j82|l . See fig. CHI By the reality property of t\ (x) and the diagrammatic symmetry 
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4k+l 



4k+l 



2k 



Figure 16: The decomposition of '7(iofc+2,6fc+2,6fc+i)/(4fc+i,4fc)(a;) into 3 parts. The hatched 
part reduces to a product of 03 (a;) 



(ED), one finds 



(10fc+2,6fc+2,6fc+l) /(4fc+l,4fc) 



X] 



^1 ='= 2L+4^) 

(4' 



X ± 



2L+4' 



k : odd 
k : even 



where the mod relation has been apphed. On the other hand, the quantum Jacobi Trudi 
formula concludes 



^10fc+2,6fc+2,6fc+l)/(4fc+l,4fc) {x) — ^(10fc+2,6fc+2,6fc+l)/(4fc+l,4fc) {x) 

which coincides with tf\x) {{tf^Y{x)) for k odd (even) in (jH^ . Thus by comparing above 
due expressions, one concludes 



tf\x± 



t)=tf\x)+t'^>ix) 



(6), 



2L + 4 



irrespective of k even or odd. 
E.6 the proof for t^^\x ± 



(4), 



2L+r 



We treat the case k even. First note the following representation of t[^^ {x + 



(112) 



1 



2L+4 



1 



9. 



t?{x + ^j^)<P{x + {6k+^-)z) 
A(4fc-i,4fc-i)/(4fc-2)(a; + 2i)T2k{x - {2k + 2)i) - k{Ak+iAk+i)/{Ak){.x)T2k-2{.x - 2ki). 
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This relation is natural by its diagrammatic representation, and can be verified easily 
by (fT^ . Again we use the reality property of t^i\x) and the diagrammatic symmetry (PT| 

to evaluate t^^\x ± ai+iO'/'l^ ^ (6^ + |)^)- 

After recombining the products of terms, it is found , 

{T2k-2{x ±2i)A(^Qk+lfik+l,4:k+l,l)/{6k,2k){x) + (T2/C (x))^A(6fc+l,6fc+l,4fc-l,l)/(6fc,2A:+2) (a^) 

—T2k{x)T2k~2{x — 2i)A(^Qk+l,6k+lAk-l,l)/(Gk,2k){x) 

—T2k{x)T2k~2{x + 2z)A(6A;+l,6fc+l,4fc+l,l)/(6A:,2fe+2)(a^)} 
+ {T2k-i{x)T2k{x) - T2k-2{X ± 2i)}T2fc(x)A(6fc+2,6fc+2,l)/(6fc+l)(a;). (113) 

Thanks to (|T9|l . the content of the first curly bracket is equal to 



T2k{x)T{i^kfik,ik,ik-l,l)/(%k-l,2k+l,2k){x) 

+Tn,,_,{x)K(L,i){x - {2k + 2)i)K(L,L)KL^i){x + {2k + 2)*). (114) 
The first term in the above is decomposed into three parts due to the tableaux rule 

2L+4' 

1 



as in the case of tf\x ± o/rriO- analogous argument leads to its expression. 



{x ± {2k _Tj^^2fc(x)A(2fc+i,2fe+i,i)/(2fc)(a; ± %{Ak + 1)) 



2' 

= (l){x ± (6A; + '^-)i)T2k{x)TD,,{^ ± (4fc + l)i) = 0(x ± {Qk + 

where Lemma El and the mod P relation are used in the second equality. 
The second term in ()114|) is identified with 



{x±t{Qk + l)t)Tr>,,_,{x)tf\x± ^ 



2' ' ' ' ^ 2L + 4: ' 

by the reality property of tf^ {x) and (j7H) . Substituting the established relation, ■* {x ± 
2lT4^) = 4^''(^) + tf\x), we obtain the following expression for the first curly bracket 
term in ()113l ) , 

0(x ±t{Qk+ ^-)^) (t^) (x) + T^,,_, (x) (tP (x) + (x)) ) 

= 0(x ± z(6fc + ^)^)(t?^(x) + tf^(x)tf^(x) + t^^\x)tf\x)). 
The second curly bracket term in pi3l ) is easily found to be 

9 

-TD2k-2i^)T2k{x)A{6k+2fik+2,i)/{6k+i){x) = -0(x ± i{Qk + -)i)TD2^_^{x)T2k{x)T^^^{x) 

= -<P{x ±i{<6k +'^)i)tf\x)tf {x). 

Lemma El is applied first, then the definitions in (f7H|) and ()89|1 is utilized. Combining 
these results, we obtain tf\x ± 21+4 = tf\x)tf\x) + t^2\^)- 
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E.7 the proof for t''^\x±^^i) = tf\x±^^i)tf\x)+tf\x)tf\x] 
Here we set k = odd. Since t^^\x) = TB.i{x)t2\x), the assertion is equivalent to 

tf{x ± ^^z) = tf{x){tf\x ± ^^.) + t^\x)Ts,{x)). 

To prove this, we prepare a few facts. 
Lemma 8. 

, / '/ \ .\ ■^{m,m~l)\X) , / \ '/ \ c; \ .x ^(m,m) I (m— 1) \ X ~\~ i) 

(f){x - {m + ^ ' ' 0(x + (m+|)z) ^ ' 

1 



3...Aim,l){x) = 5.,. A(m,m)/(l)(x + z) 



Both of them are easily proved by the 03 property. 
Lemma 9. 

Tn,,^^{x)Tn,S^ + {4k + 5)^) - Tn,,^^{x + 2z)Try,,_,{x + {4k + 3)z) 
= A(2.i)(x + (6A: + 4)z) ^ ^ 

0(x + (2A; + + (6A; + ^ ^ 

where 

^ := (r6fc+i(x + 2i)A(4fc+3,4fc+3)/(4fc+2)(a; - (2A; - 

- T6k+3{x + 2i)A(4fc+i,4fc+i)/(4fc)(a; - {2k - 8)2) 
It is easily checked that the Ihs is equal to 

TQk+l{x + 2i)T:^Qk+3fik+2fik+l)/(4k+3,4k+2){x + 2z) — TQk+3{x)'^(6k+lfikfik~l)/{4k+l,4,k){x + 4z). 

The diagram decomposition of this then leads to rhs of ()115p . 
Lemma 10. 

0(x+^2^ + ^ 

^ " ^(^3^^^^1jiy^(10'=.6^+2,6fc+l)/(2.,2.)(x) 

We use the two decompositions of T(^iok,6k+2fik+i)/(2k,2k){x). The comparison of these 
two leads to 



(f){x- {6k + ^)i) 

A(iofc,6fc+2,6fe+i)/(2fc,2fc)(a:) = , . (c^i, .3^-. T2k{x - {8k + 2)«)A(4fe_i,i)(x + {6k + 2)z) 

(pyx + yzrc + 2)^) 

-T6fc+3A(4fc+i,4fc)(x + {-2k + 2)i). 



53 



By the duality, the quantum Jacobi-Trudi formula and Lemma |H1 we find 

(j)(x + (2k + 

A(4fc~i,i) + {6k + 2)i) = ^ ^ ■^ ^(4fc+3.4fc+3)/(4fc+2) (x - {2k - 

(P^x — {6k + 

A(4k+iAk){x + {-2k + 2)i) = ^ ^ n^^^ ^(4fc+i.4fc+i)/(4fc+2)(a: - {2k - 3)i). 

Finally by noticing T2k{x — {8k + 2)i) = TQk+i{x + 2i), we prove Lemma ITHl 
The final lemma asserts a relation , 

Lemma 11. 

A{iok,6k+2fik+i)/(2k,2k){x-{6k + 5)i) = A(4kAk,2k)/i4k-i,i){x) = (f){x + {Ak + ^)i){t^^^y {x - ^) . 

The first equality is due to duality and the second comes from the definition. 
Now we prove the relation. We conveniently consider, 

{t^^\x ± — ^) - tf\x ± — 1— |_+^ 
^ 2L + 4 2L + 4 ^^i^^^+2i+4 

and show that it is equal to t^i\x + 2^+4)4'*'* + iT+i)- Fi^st we rewrite the difference as 

tf{x+j^)(^Tn^^^^{x)Tn,,{x + {Ak + 5y)-Tn,^^^^^^^ {116) 

then apply Lemma IHl With the expression of Tb^{x) in ()24|) . and the relation t^^\x) = 
t^\x)TB;i{x) , one concludes that ()llf)|l coincides with 

t^^\x+^ ■ 



3 V-^ "T 2L+4 J jj 

{x + {2k + ^)i) 



Therefore, the relation is equivalent to 

^). (117) 



D_ _+(4)^ 32 

11- 



0(x+(2fc + ^)i) ' ' 2L + 4' 

We will show this by rewriting the Ihs of pi7|l , thanks to Lemma EH 

1 



0(a: - {6k +Tj^^(ioM.Wi)/(2.,2.)(x) 
This can be further rewrite as, with aid of Lemma ITT] 



(x + (10A; + 



{x-{6k + \)i) 



'^{tf^y{x + {6k + h)i-i]-) 



which is readily shown to be equal to the rhs of ()117|1 for k odd. The case k even can be 
treated similarly. 
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E.8 the proof for tf\x ± = tfkfkf^ + t^^^' 



2L+4V — n n n ^ ^2 



We assume k even. The proof for the case k odd is similar. One substitutes tf'\x) 

tf\x ± jZ+l'^) " us^^ relation t2^''(x) = tf\x ± ai+i)- 

Then the assertion is transformed into an equivalent form, 



^■5 

2 



2L + 4 



We conveniently make a shift in 

t?ix + -7^¥i\x + —^t) = tf\x + —^i)tf\x ± —^i) 
^ ^ 2L + 4 ^ ^ ^ 2L + A ' ^ ^ 2L + A ' ^ ^ 2L + A ' 

-tf\x)tf\x + -^i) - tf\x + —^i)tf\x - -^i). (118) 
2 V ; 1 V 2L + A ' ^ ^ 2L + A ' ^ ^ 2L + A ' ^ ' 

We begin with the following lemma, 
Lemma 12. 

^ ^ 2L + 4 ^ ^ ^ 2L + 4 ^ ^ ^ 2L + 4 ^ ^ ^ 2L + 4 ^ 

-tf + — ^2)tf^(x-— (119) 
^ ^ 2L + 4 ^ ^ ^ 2L + 4 ^ ^ ^ 

2) (x + (2A; + 4)^)tf ) (x + ^^^)) ■ (120) 

Proof. Lemma El 

We first note a decomposition of t^i\x), 



(^ + ^^0 = + (6fc + 9)^)^2.(x)A(4.-i,4.-i)/(4.-2)(x + (2 A; + 4) 

-T2,_2(x + 20tf)(x-^^0 
where we have used the diagrammatic symmetry. 



1 



{x+{6k + l)i 
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Q^^{4k+i,4k+i)/(4k){x + (2A; + 2)z). 



Then the product in the Ihs of ()119|) is rewritten as, 



= 77 — , (ai I 9x-x ^2fc(a^)A(4fc-i,4fc-i)/(4fc-2)(a^ + (2/i; + A)i)tf\x + / ) 
+ (6fc + D?) 2L + 4 

-T6fc+3(x - (8/c + 2Y)tf\x + ^^^)tf)(x - (121) 

where the duahty relation T2k-2{,x + 2z) = TQk+?j{,x — {8k + 2)i) is apphed. Note that 
tf\x+ 2l+4^) = A(6fc+2,2fc)/(i)(a;)- Then we apply the decomposition rule to the second 
term in (ifeH . 

-nk+z{x-{8k + 2)i)tf\x+ ^ 



2L + 4 ' 

tf\x + 2L^i ^ ~ ^2fc(a;)A{6fc+2,2fc+2)/(i)(a;) 

, 3z (6), ^ . (6). . 2i 



—T2k (a;) A(6A,-+2,2A:+2) /(I) (a^) • 

The first equality follows from the quantum Jacobi-Trudi formula. Then the duality and 
the established relation in IE.5I lead to the last line. By substituting this into ()121|1 . one 
verifies Lemma IT^ □ 

From Lemma ^1 and ()118|) , one only has to show 

i? = tfVx + z— (x). (122) 
Our argument consists of two further steps. First we will show a lemma 
Lemma 13. 

^ ^ T2k{x) 

0(x - (2k + |)i)0(x - (6A; + |)2)0(x + (6A; + |)2) 

xA(i,i)(x + {2k + \)%) (Ti{x - {10k + 5))i)A^2k,2k)/i2k-i){x - {8k + 2)z) 
-To{x - {10k + 6)0A(2fc_i,2fc-i)/(2fe-2) {x - {8k + l)z)) . (123) 
Then the rhs of the above will be shown to be equal to that of 
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Proof. Lemma IT^ 



The definition of R in ()12n|l is not convenient in applying decomposition rules. 
We utilize the duality and rewrite 



tf\x 



A(8kAk)ix - i2k + 2)i) 



A(4fc-i,4fc-i)/(4fc-2)(a; + {2k + 4)i) = A(8kAk+2){x - {2k + 2)i). 

Then they are located at right positions to apply the quantum Jacobi Trudi formula. 
The resultant expression reads 



R 



T-2k{x) 



{x + {6k + 



X \ T8k{x - {2k + l)i)T(^8k+iAk+iAk)/(2k,2k-i){x - 2ki) 

-T8k+i{x - {2k + 2)i)T(^8k,AkAk-i)/{2k-i,2k-2){x + {-2k + 

T2k{x) 
<P{x + {<6k + 

X {Ti{x - {IQk + b)i)T(8k+iAk+iAk)/{2k,2k-i){x - 2ki) 

-To{x - (lOfc + Q)i)T,^8kAkAk-i)i{2k-i,2k-2){x + {-2k + 



where the duality is applied in the last equality. 

The "height 3" diagrams are decomposed into smaller pieces due to the tableaux rule, 
see fig. IT71 



4k+l 



'J 



2k 



2k 



Figure 17: The decomposition of the "height 3" diagram corresponding to 

'^8fc+l,4fc,4fc-l)/(2fc,2A;-l) 
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Accordingly, R is represented by 



^ _ T2fc(x)A(4fc+i,i)(x + {2k + 

0(x + (6A; + - {2k + l)i)(pix - {6k + i)z) 

X (Ti{x - {10k + 5)i)A(2fc,2fc)/(2fc-i)(a; - (8A; + 2)^) 

-Tq{x - {10k + 6)i) A(2fc-i,2fc-i)/(2fc-2) {x - {Sk + l)i) j . 
Thus the lemma is established. □ 



The final step starts from rewriting the first term in the bracket of ()123|1 in Lemma 
ITHlas. 



A{i,i)(x- {I0k + 5)i] 

Mx-{l0k + 5)i) 
_ (l){x - {10k + 
~ ^{x-{10k+f)i 



A(2fc,2fc)/(2fc-i)(a; - (8A; + 2)i) 
-^(2k+i,2k+i,i)/(2k){x - {8k + 4)i) 



+To{x - {10k + 6)i)A(2fc_i,2fc-i)/(2fc-2)(a; - {8k + l)i). 
The second term in the rhs coincides with that in the bracket of p23|) . Thus 

^ ^ T2k{x)tf\x + 

0(x - (6A; + - (lOA; + f)i) 

xA(2fc+i,2fc+i,i)/(2fc)(a; - {8k + 4)z) 

where (f){x — {10k + = (j){x + {6k + |)z) and the definition of tf^ is used. Finally 
we apply Lemma El to the above and reach the expression ()122|) . Thereby we prove the 
functional relation. 



F A result for an exceptional case 

We sketch the result for exceptional cases, L = 3 ~ 6. 
For L = 3, 4, 6, the Y system is neatly described by 



2L + 4' 



where the corresponding Dynkin diagrams are those for Er, and r = 8, 7, 6, respectively. 

The solutions of Y in terms of transfer matrices have been reported in [22] for L = 3 
and [23 for L = 4, 6. 
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The case L = 5, corresponding to M5 6+0i,2 takes a little bit complicated as announced 
in j73]. Below, we shall summarize the result then supplement outline of the proof as 
promised there. 

In this case, there are more breathers in the system than the cases we discussed above. 
We, therefore need to introduce more s/2 origin objects. 



Tb,{x) := Ti(x) (124) 

13 12 

TbM := A(8,i)(x + -O/0(x-yO (125) 

3 

TB.ix) := A(i5,8,8)/(7,7)(a;)/</)(x ± -?) (126) 

11 12 9 

Tb,{x) := A(i5,i5,8,8)/(i4,7,7)(a; + - yO0(a;± ^^)) (127) 

Tb,{x) := Tj{x) (128) 

25 

r(6)(x) := A(s,7)/(6)(x + -^) (129) 

then the following relations hold. 

3 11 12 

TB,ix±—i) = To{x ± —i) + ^{x - —i)TB.,{x) 



3 8 

TBs{x±—i) = To{x)To{x ± —i) + Tb,{x)Tb,{x) 

3 3 5 

TBsix±Y^) = To{x±j^)To{x±—i)+TB.,ix)TB,{x) 

TB,ix±^^) = To(x±i-^)+T(6)(x) 
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For a magnon-like t— system we assume, 







X — 


li' 




(130) 


''2 




'x — 


i , 




(131) 


''3 




x — 


i , 
14 




(132) 






X — 


i 




(133) 


^(2) 
''2 




x — 


i 

14 




(134) 


^(2) 
''3 




X — 


i 

u) 




(135) 






[x — 


i 

ia' 


= (x) +tS'^(x)tf^(a;)tf^(a;) 


(136) 


A3) 




X — 


i . 
14'^ 


= t^^\x)t^^\x)t^^\x) 


(137) 






'x — 


i 

14 




(138) 






x — 


i . 
14^ 




(139) 




(a) 
5i 


X ± 


i , 
14^ 


- gi''\x), a -1,2 





Remark: They are obtained from the T-system for D4, by requiring = = . 
The identification between Y and T is given by 



12 11 

Ybi = (p{x - —z)TBsix)/To{x±—i) 

Yb, = TB,(x)TB,{x)/iTo{x)To{x±P)) 

3 5 
Yb, = TBs{x)TBj{To{x±—i)To{x±—i)) 

Yb, = T^'\x)/T,{x±^^i) 



Yk,{x) = 
Y^'-\x) = 



tf\x)g?{x) 

tf\x)gf\x) 

t'^^\x)t'i\x)t^^\x) 
t^\x) 



And the explicit solution is asserted as follows. The following choice of t'f}^{x) solves 
the kink T-system 

t^^\x) = T^''\x) (140) 

tf\x) = Tb^x) (141) 

tf{x) = A(i2,8,7)/(5,4)(a;) (142) 

t^\x) = A(5,i)(x+^0Mx-^z + ^z) (143) 

t«(x) = Tb,{x)Tb,{x±'-) (144) 

4')(^) = Tb,{x±'-)Tb,{x) (145) 

tf(x) = Tb,{x ± ^)TbAx ± ^) (146) 

4')(x) = rB,(x)rB,(x) (147) 
with additional constraints, 

t^\x) = tf{x), 



^^3 - ^(3) 



g?(x) = To(x), g?{x)^To{x±j). 
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By these choice the above Y satisfy the Y— system proposed in 



Ybs 

Yk2 
Yk, 

y(2) 



, 3 . 
X ± — I 
14 

, 3 . 
14 

, 3 . 
X ± — I 
14 

, 3 . 
^=^14^ 

X ± — i 
14 

1 

X ± — i 
14 

1 

^^14^ 



Hb3(x) 

SB,(x)/:«(a;) 



We first remark a few lemmas concerning about equivalent expressions for transfer 
matrices 



Lemma 14. Tb^x) has the following two equivalent expressions in addition to M2f^) . 
Tb,{x) = Td^x) = 



The first equahty has aheady appeared in pKl ) . The second is shown by the duahty. 
Lemma 15. Tb^{x){= tf\x)) is also written as follows. 



Tb,{x) 



1 



/ 31 ., 

^^_^ipyA(4,3,I)/(2)(x--.) 



To prove this, we utihze ^i5,i5,8,8)/(i4,7,7)(3;)- The quantum Jacobi-Trudi formula and 
the duahty relation conclude, 

^(15,15,8,8)/(14,7,7)(a;) = A(i5^i5^8^8)/(14,7,7) (a^) 

+ Tq{x)Tq{x - 2i) (Ti(x + m)Ti{x -i)- To{x + 16i)T3(x)) . (148) 
On the other hand, the same quantity is given by 



'^15,15,8,8)/{14,7,7)(a:) 



0(x + '-)(P{x + ^)T7{x - llOA(8,8)/(7) {x + m) (149) 
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thanks to the semi- standard condition. Again, by the quantum Jacobi-Trudi formula and 
the duahty relation, A(8^8)/(7)(x + lOi) = Ti{x — 7i)T2{x + ^i) —Tq{x — ^i)Ti{x — 2i). Using 
this expression, as well as with some rearrangement, we find 

g 

T7(x - lli)A(8,8)/(7)(x + lOi) = Ti{x - 7i)TD^{x - i) - To{x - -i)A(2,i)(x) 

g 

+ To{x - i){T4^{x - i)Ti{x - 7i) - To{x - -i)T3{x)). 

Substitute this into ()149p , compare the resultant expression with p48|) and we have, 

i 13i 8 

A(i5,i5,8,8)/(i4,7,7)(a:) = 0(2; + -)0(a; + —){Ti{x - 7i)TD^{x - i) - To{x - -i)A(2,i)(a;)). 

By the representation of Ti(x — 7i) by 2 x 1 table, the diagrammatical rule immediately 
leads to 

A(15,15,8,8)/(14J,7)ia;j = " ^ 

where ds^ix) specifies the denominator in ()127|) . This proves the lemma after a shift in x 
and the simplification in the coefficient. 

Lemma 16. t^i\x) can be also represented by 

h (x) = ^^^A(™)/(6,2)(x + y.). 

Just as in the proof of Lemma IT31 two expressions for T(^i2,8,7)/{5,4)i^) leads to 

tf\x) = A(i2,8,7)/(5,4)(a;) = 0(3;^ 3^) (^(5.1) + §0 A{5,5)/(4) + 16i) 

-<j){x±h)Tr{x)TD,{x)). (150) 
The product A(5 + 82)A(5^5)/(4)(x + 16i) decomposes as follows. 

(A(5,i)(a; + 8i)A(5,5)/(4)(x + 16i)) = (A(5,i)(x - 4i)A(5,5)/(4)(x + 4i))^_>^+^i 

12 12 12 

= A(7,7,5,i)/(6,2)(a; + y«) + T2(X + y2)A(8,8,l)/(7)(a; + yi) 

12 3 

= A(7j,5,i)/(6,2) {X + y + T-j{x)Tb,, ± -i) (151) 

where we have used Lemma IT^ in the last equality. Using 1)1511) in ()15()j) and remembering 
Tb^{x) = Td^^x), one easily establishes the lemma. 
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Lemma 17. The following equality can be easily verified, 

3 i 12 

A(6,6,4,3,i)/(5,3,2)(a;) = 0(a; ± —i)TB,{x ± - + yz). 

The diagrammatic decomposition yields, 

A(6,6,4,3,l)/(5,3,2)(a;) = y-|^A(3,3,l)/(2)(a; ± 5i). 

The m = 3 case of Lemma El and Tb^X^) = T^^lx) then vahdates the lemma.. 

With these preparations, we prove the magnon-like t system. 
Proof of (HSO ) 

With ()14()|1 . ()129|1 and using the duality and the quantum Jacobi-Trudi formula, we can 
express the product t^i\x±Ygi) using only To(x) and T2{x). Similarly tf' can be rewritten 
by To{x),T2{x) and T4^{x), according to p42|) . Then it is immediately seen, 

To(x)tf)(x) -ti')(x± ^0 = T2{x±^t)TD,{x). 

Then the equality Tb^Xx) = T£)Xx),T2{x + yz) = T^x) and the definition of t2\x) in 
(dH) concludes (pHl l. 
Proof of (USD ) 

Consider the product A(4_3 i)/(2)(x — 4i)A(4 4_2)/(3,i)(a; + 4i). which is equal to 

<P{x±^^)XX\x±U + ^t) (152) 

due to Lemma ^1 

As is usual, the diagrammatical argument leads to 

3 5 

A(4,3,l)/{2)(a;-4i)A(4,4,2)/(3,l)(^+4«) = 0(a;±-i)0(x±-i) A(7,7,5,l)/(6,2) (a^) +r2 (x) A(6,6,4,3,l)/{5,3,2) (a^) . 

Then we apply Lemma ^1 and Lemma El to the first and the second term in the rhs, 
respectively, to reach 

0(X ± ^^i){tf\x + y2)ro(x ±h+ yZ) + T^, (x + y^Ts^^X ± ^2 + y^). (153) 



By equating ()152|) and ()153|) . we show the validity of ()131l ) . 
Proof of (HSl ) 

This is actually equivalent to ()15()j) . Note that A(5_i)(x + Si) = 0(x + ^i)t^i\x + j^i) 



and that the diagrammatic symmetry leads to A(5^5)/(4)(x — Si) = 0(x — |z)t['*''(x — ^ 
Note also T7(x) = Tb2{x). Then one easily verifies p33l ) using t2'^-'(x)in ()147p . 
Proof of (Uni ) 

This is the final nontrivial relation. This relation can be simply shown by using the 



result for the A4k+i case. We simply put k = 1 and identify t^^\ ^ tf\ t^^\ t''2\ 9i^\ 9?^ 
here with tf\X^\tf\tf\X2\t^i\tf^ there and use the duality T2 = T^. 

All remaining relations can be checked easily. Thus the above argument completes the 
proof of the magnon-like t system of the dilute A5 model . 
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G Numerical Results for the dilute Ag model in 
regime 2 

We supplement some numerical data supporting the conjecture of the analyticity of Y 
functions in section [Tl . II Throughout this appendix, we choose /3 = 0.4 and J = 1. 

The solution to the Bethe ansatz equation, corresponding to the largest eigenvalue of 
QTM, assumes the form of two strings. We list the explicit locations of roots for = 10, 

q = a 

I BAE Roots q = 

0.13548013966460354 ± 0.4464809553245548 i 

0.03263577725129405 ± 0.4935687693509952 i 

0.002550182914375612 ± 0.4961489632802383 i 

-0.02495383965703369 ± 0.4946627475810687 i 

-0.09400161546158932 ± 0.473899050451243361 
Their imaginary parts are approximately equal to ±|, irrespective of L,N,q,p. See 
figure UHl . 
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Figure 18: The Bethe ansatz roots for iV = 10 for q = 0, 0.1, 0.2, 0.3, 0.4, 0.5. Their location 
changes only slightly. 

The zeros for symmetric fusion QTM, T^, (m = 1, ■ ■ ■ , 7), are plotted in fig. (fTIHl - 
fig. ((221) • These figures clearly show that the simple-minded choice of Y in El is not 
appropriate: the choice requires should be analytic in the strip G [—1, 1]. (Or the 
arguments can be simultaneously shifted by half period so that strip is G [0.8, 1.8] U 
[—1.8, —0.8] Obviously, neither choice of strip results set of free from zeros. (Especially, 
T3 AND T5 break the analyticity in both strips.) Therefore, the transformation of Y system 
to TBA is not legitimate. 

The difference in the argument x appearing in the breather T, s/3 T and magnon-like 
t is much smaller than unity, which is crucial. 

Tb^{x) is identified with Ti{x). The functional relation then implies that it should be 
analytic within '^x G [— ^]- This can be easily checked from the figure IT9l Among other 
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Figure 19: Zeros of Ti(x)(left) and T2{x) (right) for N = 10 and q = 0.1, 0.2, 0.3, 0.4, 0.5. 
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Figure 20: Zeros of T3(a;)(left) and T^ix) (right) for iV = 10 and g = 0.1, 0.2, 0.3, 0.4, 0.5. 
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Figure 21: Zeros of 75(2:) (left) and T^ix) (right) for AT = 10 and g = 0.1, 0.2, 0.3, 0.4, 0.5. 
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Figure 22: Zeros of Tr{x) for N ^ 10 and q ^ 0.1, 0.2, 0.3, 0.4, 0.5. 
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Tm, T^ix) and T2{x + |z) are identified with t^i\x) and t'^\x). Their ANZC property in 
wider strip than '^x G [— |, |] is clearly seen. 

The breather related QTM, Tb.^{x) is equal to t^P{x), and Tb5(x)(= T£,^{x)) is 4'^'*(x), 
thus the information on their zeros is significant. Especially the latter appears frequently 
in the solution of other t^m , {m = 2, 3). 

Their zeros are shown in fig. 
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Figure 23: Zeros of Teg (x) (left) and Tbs{x){= T^^) (right) for iV = 10 and q 



0.1,0.2,0.3,0.4,0.5. 



The remaining zeros of QTM to be specified are those for t^^\x)^t'^i\x)^tf\x) and 



A (3)/ 



tf\x). They are depicted in fig. and fig.(j2Sl), respectively. 
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Figure 24: Zeros of tf^(x)(left) and t\"\x)) (right) for iV = 10 and g = 0.1, 0.2, 0.3, 0.4, 0.5. 



.(2) 



By combining the above data on zeros, one can verify that all Y and 1 + F possess no 
zeros or singularities from QTM in their own strips. Only trivial zeros (or singularity) of 
order exists for Yb^ simply coming from normalization factor 0^(x)) (or To(x)). 



H The kernel matrices for L = 5, 7 

For L = 5, we define the order as 

* {£gYB, , 1^ , ^^Yb, , hgYB, , l^gYK, , 1^ F(^) , 1^ F(') ) 
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Figure 25: Zeros of tf^(a;)(left) and tf'{x)) (right) for iV = 10 and g = 0.1, 0.2, 0.3, 0.4, 0.5. 
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For L = 7, we define the order as 



then 



68 



M = 





k 


-1 























\ 


-1 


2ch^k 





-2ch^k 





-1 






















2ch^k 


-1 

























— 1 


-1 


2ch^k 




















n 
u 







n 
u 


n 
u 


2cJi^ 


fC u 


1 

— i 


u 








n 




U 


n 


n 


n 


Z/ Uii A/ 


n 


_i 

± 








n 




(J 


n 


n 


_1 


n 




n 

u 


-1 





n 







n 


n 


n 


_1 


n 




-1 





n 







n 


n 





n 


_1 


_1 


2ch^k -1 






U 


n 


n 


n 


n 


n 


n 


-1 


2ch^^kJ 


f 

























o\ 




1 





2ch^ 


k 


1 
























1 


-2ch 


^/c 




-2ch^k 


-1 


-2ch^k 


-1 




1 


1 








-2ch 


3 7 
18 


-1 


-2ch^k 


-2ch^k 


-1 




U (J 


-1 





(J 


(J 


i 


(J 


(J 













-1 













1 



















1 













1 



















1 










1 
























1 


1 





1 




1^0 






















1 


^) 





References 

[1] G.E. Andrews, R.J. Baxter and P.J. Forrester, J. Stat. Phys.35 (1984) 193 . 
[2] S.O.Warnaar, B. Nienhuis and K. A. Seaton, Phys. Rev. Lett. 69(1992) 710. 
[3] S.O.Warnaar, B. Nienhuis and K. A. Seaton, Int. J. Nod. Phys.B 7 (1993) 3727. 
[4] A.B. Zamolodchikov, Adv. Stud. Pure. Math. 19 (1989) 641. 
[5] A.B. Zamolodchikov, Int. J. Mod. Phys. A 4 (1989) 4235. 
[6] M. Henkel and H. Saleur, J. Phys. A 22 (1989) L513. 
[7] P.G. Lauwers and V. Rittenberg, Phys. Lett. B 233 (1989) 197. 
[8] I.R. Sagdeev and A. B. Zamolodchikov, Mod. Phys. Lett. B3 (1989) 1375. 
[9] G. von Gehlen, Nucl. Phys. B330 (1990) 741. 
[10] V.P. Yurov and Al. B. Zamolodchikov, Int. J. Mod. Phys. A6 (1991) 4556. 



69 



[11] p. Dorey, "Hidden Geometrical Structures in Integrable Models" in integrable quan- 
tum field theories, NATO Advanced Study Seies, Physics V.310, eds. L. Bonuro et al 
(Plenum, New York, 1993) 



[12: 

[is; 

[14 

[is: 
[16: 
[17: 

[is: 

[19 
[20: 



[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 

[3o: 

[31 
[32 



G. Delfino and G. Musssardo, Nucl Phys. B455 (1995) 724. 

M.T. Batchelor, V. Fridkin and Y.K. Zhou, J. Phys. A, J. Phys. A29 (1996) L61. 
M.T. Batchelor and K.A. Seaton, J. Phys. A 30 (1997) L479. 
M.T. Batchelor and K.A. Seaton, Nucl. Phys. B520 (1998) 697. 
C. Korff and K.A. Seaton, Nucl.Phys. B636 (2002) 435. 

S.O.Warnaar, P.A. Pearce, B. Nienhuis and K. A. Seaton, J. Stat. Phys 74 (1994) 
469. 

B. McCoy and P. Orrick, Phys. Lett. A 230 (1997) 24. 

V.V. Bazhanov, O. Warnaar and B. Nienhuis, Phys. Lett. B 322 (1994) 198. 

U. Grimm and B. Nienhuis, "Scaling Properties of the Ising Model in a Field", 
in: Symmetry, Statistical Mechanical Models and Applications, Proceedings of the 
Seventh Nankai Workshop (Tianjin 1995), edited by M.L. Ge and F.Y. Wu, World 
Scientific, Singapore (1996), pp.384. 

U. Grimm and B. Nienhuis, Phys. Rev. E 55 (1997) 5011. 

J. Suzuki, Nucl Phys B528 (1998) 683. 

J. Suzuki, Progress in Math. 191 (2000) 217. 

P. Dorey, A. Pocklington and R. Tateo, Nucl. Phys. B661 425. 

P. Dorey, A. Pocklington and R. Tateo, Nucl. Phys. B661 464. 

F. A. Smirnov, Int. J. Mod. Phys. A6 (1991) 1407. 

L. Chim and A. Zamolodchikov, Int. J. Mod. Phys. A 7 (1992) 5317. 

A. Koubek, M. J. Martins and G. Mussardo, Nucl. Phys. B368 (1992) 591. 

A. Koubek, Int. J. Mod. Phys. A9(1994) 1909. 

P. Dorey, R. Tateo and K.E. Thompson, Nucl. Phys. B470 (1996) 317. 

Al.B. Zamolodchikov, Phys. Lett B 253 (1991) 391 

A. Kuniba and T. Nakanishi, Mod. Phys. Lett. A7 (1992) 3487. 



70 



[33] R. Tateo, Int. J. Mod. Phys. AlO (1995) 1357. 

[34] P.P. Kulisk, N. Yu. Reshetikhin and E. K. Sklyanin, Lett. Math. Phys. 5 (1981) 393. 
[35] A. Kliimper and P. A. Pearce, Physica A183(1992) 304. 

[36] A.Kuniba, T. Nakanishi and J. Suzuki, Int. J. Mod. Phys. A9 (1994) 5215, ibid 5267. 

[37] Y. Hara, M. Jimbo, H. Konno, S. Odake and J. Shiraishi, J. Math. Phys. 40 (1999) 
3791. 

[38] U. Grimm, P. A. Pearce and Y.K. Zhou, Physica A 222 (1995) 261 
[39] V.V. Bazhanov and N. Yu Reshetikhin, J.Phys. A23 (1990) 1477. 
[40] A. Kuniba and J. Suzuki, Comm. Math. Phys. 173 (1995) 225. 
[41] J. Suzuki, Phys. Lett. A195 (1994) 190. 
[42] Z. Tsuboi, J. Phys.A 30 (1997) 7975. 

[43] A.G. Izergin and V. E. Korepin, Comm. Math. Phys. 79 (1981) 303. 

[44] M. Suzuki, Phys. Rev.B.31 (1985) 2957. 

[45] T. Koma, Prog. Theoret. Phys. 78(1987) 1213. 

[46] M. Suzuki and M. Inoue, Prog. Theor. Phys. 78 (1987) 787. 

[47] J. Suzuki, Y. Akutsu and M. Wadati, J. Phys. Soc. Japan 59 (1990) 2667. 

[48] M. Takahashi, Phys. Rev. B 43 (1991) 5788, see also vol. 44 p. 12382. 

[49] A. Kliimper, Ann. Physik 1 (1992) 540. 

[50] A. Kliimper, Z. Phys. B 91 (1993) 507. 

[51] G. Jiittner, A. Kliimper and J. Suzuki, Nucl. Phys. B512(1998) 581. 
[52] A. Kuniba, K. Sakai and J. Suzuki, Nucl Phys B525 597. 
[53] K. Sakai and Z. Tsuboi, Int. J. Mod. Phys. A15 (2000) 2329. 
[54] T.R. Klassen and E. Melzer, Nucl. Phys. B350(1991) 635. 
[55] Y.K. Zhou, unpublished note 

[56] I. Krichever, O. Lipan, P. Wiegmann and A. Zabrodin, Comm. Math. Phys. 188 
(1997) 267 

[57] R.M. EUem and V.V. Bazhanov Nucl.Phys. B647 (2002) 404. 

71 



[58] R. Kedem, T.R. Klassen, B.M. McCoy and E. Melzer, Phys.Lett. B307 (1993) 68-76. 

[59] S. Dasmahapatra, R. Kedem, T.R. Klassen, B.M. McCoy and E. Melzer, 
Int.J.Mod.Phys. B7 (1993) 3617. 

[60] A. Berkovich, Nucl.Phys. B431 (1994) 315. 

[61] E. Melzer, Lett.Math.Phys. 31 (1994) 233. 

[62] A. Berkovich and B.M. McCoy, Lett.Math.Phys. 37 (1996) 49 . 

[63] S. O. Warnaar, J. Stat. Phys.82 (1996) 657. 

[64] S. O. Warnaar, J. Stat. Phys.84 (1996) 49. 

[65] S.V. Kerov, A. N. Kirillov and N. Y. Reshetikhin, J. Sov. Math. 41(1988) 916. 

[66] A. N. Kirillov and N. Y. Reshetikhin, J. Sov. Math. 41(1988) 925. 

[67] G. Hatayama, A. Kuniba, M. Okado, T. Takagi and Y. Yamada, Contemporary Math. 
248 (1999) 243. 

[68] G. Hatayama, A. Kuniba, M. Okado, T. Takagi and Z.Tsuboi, Prog. Math. Phys. 23 
(2002) 205. 



[69] A. Schilling, "Rigged configurations and the Bethe Ansatz" (math-ph/0210014) 



[70] A. Berkovich, B.M. McCoy and P.A. Pearce , Nucl.Phys. B519 (1998) 597 . 

[71] A. Berkovich and B.M. McCoy, "The perturbation 02, i of the M{p,p + 1) 
models of conformal field theory and related polynomial character identities" 
(math.QA/9809066|). 

[72] S.O. Warnaar and P.A. Pearce, J.Phys. A27 (1994) L891. 

[73] J. Suzuki, J. Phys. 37 (2004) 511. 



72 



